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ABSTRACT 


This  study  deals  with  several  two-dimensional  scattering  and  diffraction 
problems  in  anisotropic  media.  The  intent  is  twofold:  First,  to  generalize 
mathematical  methods  applicable  in  isotropic  regions  to  a  certain  class  of  aniso¬ 
tropic  problems;  and  Second,  to  study  the  solutions  of  the  anisotropic  problems  in 
3uch  a  manner  as  to  highlight  certain  common  properties  which  point  the  way  toward 
the  construction  of  approximate  solutions  for  configurations  with  more  general 
structural  shape  or  anisotropy.  In  the  low-frequency  range,  the  method  of  multipole 
expansion  is  used.  It  is  demonstrated  that  the  rigorous  solution  for  the  problem 
cf  scattering  of  a  plane  or  cylindrical  wavd  by  an  obstacle  which  is  small  compared 
to  the  wavelength  may  be  expanded  in  a  series  whose  terms  correspond  to  multipole 
radiation  in  the  anisotropic  medium.  As  an  illustration,  the  excitation  coefficients 
of  the  first  few  terms  arising  from  scattering  by  a  narrow  conducting  ribbon  are 
calculated.  In  the  high  frequency  range,  geometrical  optics,  the  first-order 
asymptotic  solution  of  Maxwell's  equations  is  considered  first.  The  ray  refractive 
index  is  calculated,  and  the  laws  of  propagation  and  reflection  of  rays,  which  define 
the  trajectories  of  energy  transport,  are  derived  from  Fermat's  principle.  To 
obtain  an  insight  into  diffraction  phenomena,  two  types  of  representative  problems  -- 
--  diffraction  by  a  straight  edge,  and  diifraction  by  a  smoothly  curved  object  --  are 
discussed.  Rigorous  solutions  are  found  which  are  then  expanded  asymptotically 
for  high  frequencies.  The  diffracted  field  contributions  are  phrased  in  a  manner 
which  emphasizes  the  local  character  of  the  diffraction  process  by  exhibiting  an 
explicit  dependence  on  the  local  properties  of  the  scattering  object  and  on  the  adjacent 
medium.  The  asymptotic  expressions  may  be  put  in  a  form  which  constitutes  a 
generalization  of  Keller's  geometrical  theory  of  diffraction,  thereby  providing  a 
method  for  constructing  high  frequency  solutions  for  a  rather  general  class  of 
diffraction  problems  in  anisotropic  regions. 

Special  attention  is  paid  to  the  case  for  which  the  medium  parameters  are 
complex  or  negative  real,  so  that  the  results  apply  to  a  certain  class  of  plasmas 
with  losses  or  with  "hyperbolic"  refractive  index  characteristics. 
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CHAPTER  I 

INTRODUCTION  AND  SUMMARY 


To  gain  an  understanding  of  the  phenomena  of  scattering  and  diffraction 
of  electromagnetic  waves  by  objects  embedded  in  anisotropic  dielectric  media 
has  been  the  goal  of  many  investigators  in  recent  years.  The  general  problem 
is  rather  complex,  but  various  special  problems  have  been  solved,  both  for 
their  own  sake  and  for  the  sake  of  extrapolating  thu  obtained  results  to  other, 
more  difficult  configurations.  In  the  following  investigation,  the  same  course 
has  been  chosen.  Several  special  problems  have  been  treated,  and  an  attempt 
has  been  made  to  find  the  limits  of  validity  of  the  methods  used. 

Throughout  this  work  we  assume  that  the  surrounding  medium  is  aniso  ¬ 
tropic  and  can  be  characterized  by  a  dielectric  tensor  C  the  elements  of 
which  are  known  functions  of  frequency  and  other  parameters^.  Thus,  at  a 
fixed  frequency  (steady  state),  the  tensor  elements  are  regarded  as  constant, 
independent  of  the  fields  under  investigation.  The  limitations  of  this  approach 
have  to  be  kept  in  mind.  However,  results  obt-.ined  from  such  an  idealized 
model  for  the  medium  are  known  to  furnish  workable  approximations  in  such 
studies  as  communicc.tion  through  the  ionosphere,  radio  wave  propagation 
around  the  earth  etc.  .  where  the  propagating  fields  are  not  very  intense.  Of  the 
mathematical  methods  used  in  this  work,  two  should  be  specifically  mentioned: 

L.  B.  Felsen^  devised  a  method  of  reducing  boundary  value  problems  which 
arise  from  a  certain  class  of  electromagnetic  wave  problems  in  uniaxially 
anisotropic  media  to  other  boundary  value  problems,  corresponding  to  wave 
propagation  in  isotropic  media.  This  is  accomplished  by  means  of  a  simple 
transformation  of  variables.  The  transformed  problems  may  then  be  treated  by 
known  methoas,  and  the  solutions  transformed  back  to  the  original  variables. 
Although  this  method  is  applicable  only  for  a  rather  restricted  class  of  problems, 
its  value  lies  in  enabling  one  to  find  rigorous  solutions.  These  solutions  may 
then  be  evaluated  approximately  by  asymptotic  methods  in  the  various  wave- 
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length  ranges.  The  asymptotic  results  may  be  interpreted  in  physical  terms, 
and  their  forms  indicate  how  solutions  of  more  complicated  problems  (to  which 
it  might  be  impossible  to  find  rigorous  solutions)  should  appear.  This,  in 
effect,  is  the  main  goal  of  this  investigation. 

In  finding  asymptotic  solutions  in  the  short  wavelength  ("quasi  optic") 

(3) 

range,  use  has  been  made  of  the  method  formulated  by  J.  B.  Keller  for  iso¬ 
tropic  electromagnetic  problems.  With  the  help  of  his  technique,  one  may 
construct  asymptotic  solutions  to  a  rather  broad  class  of  wave  problems.  It 
is  important  to  n^te  that  this  theory  has  not  yet  been  shown  to  hold  generally. 
Instead,  many  specific  problems  which  can  be  solved  rigorously  have  been 
shown  to  agree  with  this  theory.  In  order  to  extend  the  applicability  of  Keller's 
theory  to  anisotropic  media,  it  is  therefore  necessary  to  find  rigorous  solutions 
of  various  prototype  problems,  (with  the  help  of  Felsen's  method,  or  otherwise), 
and  to  deduce  from  their  asymptotic  behavior  the  extension  of  Keller's  method 
to  a  class  of  electromagnetic  wave  problems  in  anisotropic  media.  It  has  been 
mentioned  that  the  elements  of  the  dielectric  tensor  s  are  regarded  as  con¬ 
stants  at  any  given  frequency.  These  constants  may  be  complex  numbers  and, 
for  example,  in  a  magnetoplasma,  may  lie  on  the  negative  read  axis  of  the  com¬ 
plex  plane.  Under  these  conditions,  the  wave  equation  for  some  of  the  field 
components  becomes  a  hyperbolic  rather  than  an  elliptic  equation.  Although 
Felsen's  method  may  be  applied  formally  even  if  the  elements  of  c  are  not 
real  and  positive,  the  solutions  have  to  be  carefully  examined  to  find  whether 
they  are  still  physically  meaningful.  The  solutions  of  several  problems  have 
been  investigated  as  functions  of  the  complex  elements  of  e  and  the  regions  of 
their  validity,  has  been  determined. 

In  Chapter  II  the  class  of  physical  problems  which  is  considered  in  this 
work  is  described,  the  corresponding  boundary  value  problem  formulated,  and 
Felsen's  method  of  solution  is  outlined.  In  Chapter  III,  this  method  is  used  for 
finding  an  exact  solution  to  the  problem  of  diffraction  by  an  elliptic  cylinder. 

This  solution  is  then  expanded  asymptotically  in  the  long  wave  range  ("multipole 
expansion"),  and  the  special  case  of  a  narrow  conducting  ribbon  is  worked  out  in 
some  detail.  Chapter  IV  covers  several  aspects  of  the  short  wave  diffraction 
problem.  First,  in  order  to  utilize  ray  concepts,  geometrical  optics  is  dis- 
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cussed  and  some  general  properties  of  the  geometric  optical  field  are  derived 
for  our  class  of  problems.  Next,  diffraction  effects  are  investigated,  from 
the  point  of  view  of  Keller's  geometrical  theory  of  diffraction.  A  representative 
problem  for  the  treatment  of  edge  diffracted  rays  is  the  diffraction  by  a  con¬ 
ducting  half-plane,  which  is  solved  rigorously  by  means  of  the  Wiener-Hopf 
(4) 

technique  .  The  asymptotic  expansion  of  the  rigorous  solution  is  interpreted 
in  terms  of  geometric  optical  and  diffracted  rays.  For  the  treatment  of  surface 
diffracted  rays,  two  problems  are  solved:  diffraction  by  a  parabolic  and  by  an 
elliptic  cylinder.  Rigorous  solutions  are  found  by  Felsen'a  method,  and  again 
the  asymptotic  expansions  interpreted  in  ray-optical  terms.  The  diffraction 
coefficients  and  decay  exponents  found  in  that  way  may  be  used  to  construct 
solutions  to  more  general  diffraction  problems  of  the  same  class,  exact  solu¬ 
tions  of  which  are  unavailable.  It  is  demonstrated  that  these  diffraction  co¬ 
efficients  and  decay  exponents  are  functions  of  the  ray  refractive  index  in  the 
medium.  A  ray  refractive  index  may  be  found  for  more  general  anisotropic 
media  than  those  considered  in  this  investigation.  Thus,  the  formulation  used 
in  this  study  indicates  the  possibility  of  applying  Keller's  theory  to  more  general 
problems,  of  which  those  considered  here  are  special  examples. 
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CHAPTER  II 

FORMULATION  OF  THE  PROBLEM 


We  consider  two-dimensional  electromagnetic  problems,  in  which  the 
fields  are  generated  by  z- directed  line  currents  of  constant  strength  and  the 
scattering  obstacles  are  perfectly  conducting  cylindrical  objects  whose  axis  is 
parallel  to  z  (Fig.  1) 


Fig.  i  The  geometry  of  the  two-dimensional 
diffraction  problem. 
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The  surrounding  anisotropic  medium  is  assumed  to  be  described  by  a  dielectric 
tensor  having  the  form 


r  Z_  Z  G 
— 0 — O  z 

(2- la) 

[a  ie\ 

(2-lb) 

W8  £  / 

Initially  a,  8  and  e  are  assumed  to  be  real  and  positive  constants.  Later  on, 

analytic  continuation  into  the  complex  plane  will  be  investigated  for  some  of 

these  elements.  In  this  class  of  problems.  Maxwell's  equations  have  E-mode 

solutions  (H  =  H  =0)  which  are  excited  by  magnetic  line  currents,  and 
x  y 

H-mode  solutions  (E  =  E  =0)  which  are  excited  by  electric  line  currents. 

X  y  (14) 

The  fields  in  the  latter  category  behave  essentially  as  in  an  isotropic  medium 

and  will  not  be  considered  further. 

-ixt  t  * 


For  an  implied  e 


time  dependence,  the  pertinent  field  equations 


are 


V  x  H  =  -ix  s  •  E  (2-Za) 

V  x  E  =  ixn  H  -  M  (2-2b) 

—  o —  — 

With  the  excitation  in  the  form  of  a  magnetic  line  current  of  unit  strength 

M  =  6(x  -  x')  6(y  -  y  ')  (2-3a) 

the  magnetic  field  has  the  simple  form 

H  =zQH(x,y)  (2- 3b) 


where  H(x,  y)  satisfies  the  wave  equation 

|u,.L_  +  s£__  +  k  2  & ]  H(x,  y)  = -ixe  A  6(x-x')  6  (y-y ') 
\  dx  d}'  ° 


(2-4) 
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with 


(2- 4a) 


A 


a  e  - 


a  iB 
-i8  e 


(2- 4b) 


The  field  components 
(2-2a): 


and  E 

y 


can  be  derived  from 


H(x,  y) 


through  Eq. 


E 

x 


3H 

*7 


+  iB 


E 

y 


dH 

*7 


+  a 


A 


A 


(2-5a) 


(2-5b) 


The  boundary  conditions  to  be  satisfied  by  H(x,  y)  are  a  radiation  condition* 
at  r-*°°  and  the  vanishing  of  the  tangential  component  of  the  electric  field 
on  the  surface  of  the  perfectly  conducting  scatterer.  If  the  object  has  sharp 
edges  as  in  Fig.  1  the  field  must  also  satisfy  an  edge  condition. 


The  dielectric  tensor  defined  by  Eqs.  (2- la,  b)  represents  two  cases  of 
particular  interest:  when  a  =  e  in  Eq. (2-lb),  the  tensor  represents  a  gyrotropic 
medium.  It  will  be  shown  in  Chapter  IV  that  the  ray  refractive  index  for  such  a 
medium  does  not  depend  on  direction  (for  the  two-dimensional  class  of  problems 
considered  here).  The  only  difference  between  this  case  and  the  isotropic  case 
is  in  the  boundary  conditions  at  an  interface,  as  can  be  seen  from  Eqs.  (2-5a,  b). 


*  Since  the  medium  is  anisotropic,  the  radiation  condition  requiring  the  out¬ 
ward  flow  of  energy  cannot  be  phrased  simply  as  the  "outgoing  wave"  condition 
familiar  from  isotropic  problems.  This  aspect  has  received  detailed  attention 

elsewhere'*). 
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When  8  =0  in  Eq.  (2-lb),  the  tensor  represents  a  uniaxial  medium.  Ionized 
gases  subject  to  a  very  strong  constant  magnetic  field  in  the  y-direction  may 
be  represented  (approximately)  by  such  a  tensor.  The  frequency  dependence  of 
e  in  Eq.  (2-lb)  for  such  media  indicates  the  possibility  that  e  becomes  negative 
in  a  certain  range  of  frequencies.  If  losses  due  to  collisions  and  other  effects 
are  considered,  £  will  be  a  complex  number  with  I  m  e_^  0  for  the  implied 
time  dependence.  *  Further  discussion  will  be  limited  mostly  to  this  latter 
case  (8  =  0)  as  it  is  the  simplest  case  which  displays  the  anisotropy  both  in  the 
ray  refractive  index  and  the  boundary  conditions. 

It  is  the  object  of  this  study  to  find  solutions  of  Bj.  (2-4)  subject  to  the 
said  boundary  conditions  for  a  variety  of  configurations,  to  evaluate  asymptotic 
expansions  of  the  solutions  at  low  and  high  frequencies,  and  to  interpret  these 
in  physical  terms.  Further,  the  dependence  of  the  solutions  on  the  parameter  € 
(as  it  is  allowed  to  take  on  complex  values)  will  be  investigated  to  ?.ssure  that 
they  remain  valid  throughout  the  region  of  interest,  namely  lme>0. 

Felsen's  method  of  solving  Eq.  (2-4)  in  the  case  8  -  0  is  as  follows. 
Define  a  change  of  variables 


which  transforms  Eq.  (2-4)  into 


H(u, v)  =  -i-x£o4/aF  6(u-u')  6(v-v;) 


where 

k  =  k 

o’ 


(2-6) 


(2-7) 


(2-7a) 


*  In  reference  (1)  the  reader  may  find  a  detailed  discussion  of  the  properties  of 
the  dielectric  tensor  of  ionized  gases  under  various  conditions. 
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Fig.  2  (a)  The  real  space  (x,  y). 

(b)  The  transformed  space  (u,  v). 


On  the  surface  of  a  perfect  conductor,  E  =  0,  which  can  be  expressed  in 
terms  of  H  via  Eqs.  (2-5) 


1  n  x  E I  =  E  sin 9  -  E  cos  6  =  . — — 
—  —  y  x  i  a’  e 


dH  sin  8  +  3H  cos  0 
dx  e  a 


where  n  is  a  unit  vector  normal  to  S.  One  notes  the  relations 


=  0  onS 

(2-8) 
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dH 
Tx  - 

1  3H 
Ja  TS 

(2-9a) 

dH 

*7  = 

1  dH 
Tv 

(2-9b) 

sin  6 

=  = 

Js  dv 

ft  sin  9 

(2-9c) 

ds 

V a  du2  +  edv^ 

k  o  A 

y  a  cos^  0  +  e 

sin^  9 

cos  0 

dx 

Jo.  du 

A 

Ja  cos  8 

(2-9d) 

~  ds 

du2  +  edv2 

cos2  S  +  e 

sin2  &  * 

with  the  quantities  sin0  and  cos  9  defined  as 


sin  6  =  ■■  j-v-,  (2-9e) 

^du2  +  dv2 


cos  8  =  -========■  (2-9f) 

«ydu2  +  dv2 


Substitution  of  Eqs.  (2-9  a-d)  into  Eq.  (2-8)  yields  the  condition. 


dH 

Tu 


.  *  ,  dH  $  dH 

sin  53  +  •  COS  C  =  g 


=  0 


on  S 


(2-10) 


S  is  given  by  a  function  f(u,  v,  e,(i)  =  0.  It  is  obtained  from  the  expression  of 

A 

S  by  means  of  Eq.  (2-6).  If  a,  e  are  not  real  and  positive  S  does  not  have  the 

geometrical  meaning  of  a  surface  in  the  u-v  space.  It  is  seen  that  the  oblique 

derivative  boundary  condition  in  Eq.  (2-8)  becomes  a  conventional  Neumann 

type  condition  in  the  u-v  space  (Eq.  (2-10) ).  If  arbitrary  complex  values  of 

a  and  e  are  allowed,  the  cerivation  of  Eqs.  (2-7)  from  (2-4)  and  Eq.  (2-10) 

from  (2-8)  does  not  change.  In  this  case,  S  cannot  be  drawn  as  in  Fig.  2  (b) 

* 

because  u  and  v  are  complex.  Also  5  as  defined  by  Eqs.  (2-9  e,  f)  becomes 
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complex.  Nevertheless,  if  a  solution  is  found  which  satisfies  3q.  (2-7)  and 

A 

the  boundary  condition  (2-10)  on  S  (which  is  given  now  by  a  function  of  the  com¬ 
plex  variables  u,  v),  it  will  be  a  solution  of  Eq.  (2-4)  with  the  boundary  condition 
(2-8)  on  S  (which  is  a  real  surface)  in  the  x,  y  space.  If  such  a  solution  is  in 
the  form  of  an  integral  or  an  infinite  series,  with  a,  e  appearing  as  parameters, 
it  will  be  meaningful  only  in  those  regions  of  a.  and  e  where  it  is  convergent 
and  satisfies  the  radiation  condition.  The  radiation  condition  requires  radial 
outflow  of  energy  from  the  source  at  distances  very  large  compared  to  the 
wavelength.  In  an  anisotropic  medium,  the  wavelength  is  a  function  of  direction, 
and  the  expression  kr  =  2n  »1  should  be  replaced  by  ko  N(0)  r»  1,  where 

N(6)  =  =^xcos20  +  esin^o  (2-11) 

is  defined  as  the  ray  refractive  index^^^’,  and  M3)  is  the  wavelength  along 
a  radial  line  from  the  source,  which  makes  an  angle  8  with  the  y  axis. 

Proof  of  the  right  hand  part  of  eq.  (2-11)  will  be  given  in  Chapter  IV,  section  B . 

Using  Eq.  (2-11)  it  is  easy  to  show  that  kQN(S)  r»  1  expressed  in  the 
u-v  coordinates  becomes  kVu+v^ »  1,  which  is  just  the  usual  requirement 
in  an  isotropic  medium.  Thus  the  transformed  problem  comprises  Eq.  (2-7), 

A 

subject  to  a  Neumann  type  boundary  condition  on  S  and  a  radiation  condition 
of  the  usual  type.  This  is  a  conventional  diffraction  problem  in  an  isotropic 
medium  characterized  by  a  wavenumber  k.  If  its  solution  can  be  found,  one 
•may  use  the  inverse  transformation  of  Eq.  (2-6)  and  thus  obtain  a  solution  for 
the  original  problem  in  the  x-y  space.  As  long  as  a.  and  e  are  real  and 
positive,  the  physical  nature  of  the  origin  _  and  transformed  configurations  is 
retained,  and  a  solution  which  is  physical  in  the  u-v  space  will  be  physical  in 
the  x-y  space  as  well,  because  Eq.  (2-6)  introduces  a  change  of  scale  only. 

If  however  these  parameters  are  complex  or  if  one  is  negative  real,  it  has  to  be 
shown  that  solutions  obtained  in  a  formal  manner  via  Eq.  (2-6)  still  solve  the 
original  physical  problem,  by  checking  their  convergence  and  the  satisfying  of 
the  radiation  condition.  In  the  absence  of  physical  boundaries  (radiation  in 
infinite  space)  or  in  the  presence  of  a  perfectly  conducting  plane  boundary,  it 
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has  been  shown  that  Felsen's  method  yields  physically  meaningful  solutions 

(2) 

if  a,  e  are  restricted  to  certain  regions  in  the  complex  plane  . 
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CHAPTER  IH 

SCATTERING  IN  THE  LONG  WAVE  RANGE 


A.  FORMULATION  OF  THE  PROBLEM  OF  SCATTERING  BY  AN 
ELLIPTIC  CYLINDER 

In  this  chapter  we  solve  the  problem  of  scattering  by  a  perfectly  con¬ 
ducting  elliptic  cylinder  embedded  in  a  uniaxial  medium.  The  rigorous  expres¬ 
sion  for  the  Green's  function  is  found  by  utilizing  Felsen' s  method.  An 
asymptotic  expansion  is  obtained  which  is  useful  when  the  dimensions  of  the 
cylinder  are  small  compared  to  \o.  For  the  special  case  of  a  narrow  ribbon, 
the  terms  in  the  asymptotic  solution  are  identified  as  multipole  radiations  in 
the  uniaxial  medium. 


Fig.  3  The  geometry  of  the  elliptic  cylinder. 
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It  has  been  shown  in  Chapter  II  that  all  the  field  components  are  derivable  from 
a  Green's  function  G  which  satisfies  the  equation 


+ 


G  = 


-  6  (u-u')  6  (v-v ') 


subject  to  the  conditions 


(3-1) 


3G  _  n  e 
— *-  =  0  on  S, 
an 

T'c  2 

Radiation  condition  at  y  u  +  v  —  • 

The  equation  for  the  given  surface  S  in.  the  anisotropic  x-y  space  is 


(3-la) 

(3-lb) 


(3-2a) 


which  may  be  expressed  in  terms  of  u,  v  by  transforming  first  to  the  x,  y 
coordinates,  and  then  using  Eq.  (2-6).  The  resulting  equation  for  S  in  the 
u-v  space  is: 

(3-2b) 

(3-2c) 

(3-2d) 


A(S  )u2+B(3  )v24  2C(3  )  uv 
'  o  o  o 


=  1 


where 


A(3 


/  .  2, 

/sin  3 

j-A-r 


2,  Zn 

sin  3  cos  c 

I 


/cos  6  sin  9 

B<So>  *  5  (  T"2 


2-  .  20  \ 
3  sin  9  \ 

•-tf 


(3-2e) 


I 
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In  order  to  have  the  equation  for  S  in  a  convenient  form,  one  may  perform  a 
rotation  through  an  angle  6  in  the  u-v  space 


u  1  / cos  6  -  sin  6\  /U 
v  /  tsin  6  cos  6/  \  V 


(3-3) 


so  th  t  the  equation  for  S  in  the  U-V  space  may  be  written  as 

|A(0&)  cos26 +  B(0q)  sin26  +  2C(0q)  sin 6  cos  U2  + 

+  [A(6  )  sin26+B(0Q)  cos26 -20(0^)  sin  &  cos  &J  V^  + 

4- 2 1  |b( 0 -  A(0o)]  sin  6  cos  6  +  C(0q)  (cos26 -  sin26)  |  UV  =  1  (3-4) 

If  one  sets  the  coefficient  of  UV  equal  to  zero,  the  angle  6  is  determined  in 
terms  of  quantities  a,  b,  0  q,  a,  e , 


tan  26  = 


Zjae  sin  6  cos  6  i-4~  * 

o  o  \_2 

_ _ i  ....  — 

■  2r 
sin  S 


To  s\  .  2 0  fa  FT 

f — sr-  -  — +  cos  o  1 — s-  — 

1  2  ,21  o\.  2  2/ 

\a  by  \d  a  / 


(3-5) 


and  the  equation  for  S  now  has  the  canonical  form  of  an  ellipse  in  the  U-V 
frame; 


U2  V2 

2  +  T  "  1 

m  n 


(3-6) 


With  the  help  of  Eq.  (3-5),  m  and  n  may  be  expressed  in  terms  of  A,  B  and  C 
which  are  defined  by  Eq.  (3-2c,  d,  e).  The  results  are 


m  = 


n  = 


j  [(A  +  B)  +  4  (A  -  B)2  +  4CJ 
j  [(A  +  B)  -  /(A  -  B)2  +  4C2  ]  | 


-  1/2 
-  1/2 


(3-7a) 

(3-7b) 


Some  algebraic  manipulation  shows  that  for  real  and  positive  a  and  g  ,  m 
and  n  (the  semi -axes  of  the  ellipse  in  the  U  -  V  space)  are  real  and  positive 
quantities.  For  complex  or  negative  real  values  of  a  and  (or)  e  ,  U,  V  become 
complex  variables,  and  Eq.  (3-6)  is  no  longer  the  equation  of  a  real  ellipse. 

The  boundary  value  problem  is  now 


G  =  -  6  (U  -  u')  6  ( V  -  V7) 


(3-8a) 


Radiation  condition  at 


00 


(3-8b) 


(3-8c) 


B.  SOLUTION  IN  ELLIPTIC  CYLINDER  COORDINATES 

The  boundary  condition  (3 -8b)  suggests  solution  of  the  problem  in 
elliptical  cylinder  coordinates. 


♦  V 


Fig.  4  Elliptical  cylinder  coordinates. 
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2  ? 

Sem(!VV’0,?vtheir  linear  combination  me  (t];1i  ).  The  "radial"  Mathieu  func- 
•  (10)  *** 
tions  Mm(§;h)  are  so  defined  as  to  behave  like  outgoing  or  incoming  waves 

(or  their  linear  combinations)  as  §-*».  In  terms  of  these  functions  the  solution 

to  Eqs.  (3-8)  is  „ 

G(§,r);S  V)  =-|- Y  ^  mem(ri')  mem(-ri)  M^(%  )  - 


IM(1)  (§<) 


M7 

m 


U) 


«b> 


m 


where 


M7 

m 


UT 


m  {3)  (?  <)| 


m 


=  G .  +  G 

d  s 


§ 


> 


§  when  §  >  57 
§  7  when  §7  >  § 


(3-11) 


(3 -11a) 


§ 


< 


and 


§  when  5  <  §' 
?  'when  §7  <  § 


M 


(1) 
/( 3) 


m 


<^o> 


(1) 

(§  ;  h) 

— 


the  derivation  of  Eq.  (3-111  is  given  in  appendix  A. 
The  first  part  of  Eq.  (3-11) 
oo 

Gd  *  ?  I  mem  <"'> 
m=-oo 


(3 -lib) 


(3 -lie) 


(3-12) 


is  a  representation  of  the  primary  field,  for  which  a  closed  form  expression 
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is  known^ 

Gd  =T  Ho}  (k|  R- R'l)  H^kjiu.u^  +  (V-V^]  = 

=T  H(0U  i_  HJ)1,[koN(e)  Ir  -  r'|]  (3-13) 

A 

with  R,  R'.fp  and  cp  '  defined  in  Figure  4,  jr.jr'  and  9  defined  in  Figure  3  and 

N(§)  defined  by  £q.  (2-11).  Eq.  (3-13)  may  of  course  be  derived  directly  by 

applying  Felsen's  method  to  Eq.  (2-4)  in  an  unbounded  space.  This  will  be  done 

in  Chapter  IV,  section  B  .  The  second  part  G  of  Eq.  (3-11)  is  a  representa- 

(9-12)  8 

tion  of  the  scattered  field.  It  is  known  that  this  representation  may  be 

expanded  as  a  power  series  in  h  .  If  h«l  and  §  «1,  i.  e.  when,  the  dimen¬ 
sions  of  the  scattering  obstacle  are  small  compared  to  wavelength,  this  series 

is  rapidly  convergent^  *  The  form  of  the  scattered  field  G  becomes 

s 

quite  simple  in  the  "far  field"  region,  i.  e.  when  both  source  and  observation 
points  are  located  at  large  distances  (compared  to  wavelength)  from  the  obstacle. 
Then  h«l  and  cosh?»l,  in  which  case^*^ 

(§)~J  0<R)  (3-14a) 

m  m 

(?)~H^}(kR)  (3- 14b) 

m  m 


Because  of  the  rapid  convergence  of  the  series  expansion  of  G  ,  only  a  finite 

8  (1) 

number  of  terms  contribute  effectively,  and  it  is  permissible  to  replace  H  •  (kR) 

/  t  \  •  FT  /  O  Ki 

by  its  large  argument  asymptotic  approximation  H'  '(kR)  e  .  Thus  G  is 

O  0 

written  as 


°s' [tho  ,(kR>]  EiHo ’<kRi  = 

=  jTHo  ’  ^oN(e)r||7Ho  ’  [koN<9'>r']|i<e’e'>  <3-15) 


#  There  are  differences  between  the  definition  and  notation  of  Mathieu  functions 
in  references  (8)  and  (11).  In  this  study,  the  notation  of  reference  (8)  is  used. 
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i 


x 


where  it  is  understood  that  the  large  argument  approximations  of  the  Hank«l 
functions  are  implied  (i.  e. ,  the  result  is  good  to  G(l//kR)  and  0(  1  /^icR  ) ) , and 
where 


?  M  (1)(§  ) 

gfa.cp')  =  -4i  )  (-l)mme  (ti')me  (-»))  - — — 

m=^»  m  m  M  (3)'§  ) 


m 


(3-15a) 


g  may  be  expanded  as  a  power  series  in  h  which,  due  to  the  condition  h«l, 
will  be  rapidly  convergent.  According  to  Eq.  (3-9g),  11  and  T| '  in  this  expan¬ 
sion  may  be  replaced  by  <p  and  cp'.  All  pertinent  formulas  for  this  expansion 
are  given  in  references  (9)  an.-1  -10),  and  the  calculation  has  been  carried  out  in 
detail  in  reference  (12).  The  scattered  field,  given  by  Eq.  (3-15),  has  the  form 
of  a  cylindrical  wave  radiated  from  the  obstacle,  with  an  amplitude  proportional 
to  the  incident  field  at  the  location  of  the  obstacle  and  to  a  "scattering  coefficient" 
g(cp,cp'),  which  is  a  function  of  the  geometry  of  the  scatterer  and  the  properties 
of  the  medium.  The  identification  of  g(<p , cp  ' )  as  a  superposition  of  multipole 
radiations  in  the  given  medium  will  be  carried  out  in  the  next  section  for  the 
special  case  of  5q=0,  i.  e.  where  the  elliptic  cylinder  reduces  to  a  atrip. 

C.  SCATTERING  BY  A  NARROW  STRIP 


If  the  ellipse  in  Fig.  3  reduces  to  a  strip,  then  b=0,  and  a=cQ  is  its 
half-width.  To  find  the  half-width  m  of  the  transformed  strip  (n*0  in  this 
case),  one  may  use  Eqs.  (3-2)  and  (3-7).  The  calculation  is  rather  involved 
since  one  has  to  determine  lim  m(a,  b)  in  Eq.  (3 -7a)  in  order  to  obtain  the 

b-0 

desired  result.  A  more  straight-forward  approach  may  be  used  in  this  special 

case.  It  consists  of  transforming  the  end  points  of  the  strip,  which  are 

P,  5  (c  sin  8  ,  c  cos  9  )  and  P,  =  (-c  sin  0  ,  -c  cos  8  )  in  the  xy  space, 

lOOOO  CO  O  O  O  '  * 

into  the  uv  space  via  Eq.  (2-6).  The  result  will  be 


(3-16) 


{ 
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where  N(9q)  is  defined  by  Eq.  (2-11).  Eqs.  (3-lGa)  and(3-9a,b)  then  yield 


k  a 

h  =  -§-  N(9o), 

(3-17a) 

?o=0’ 

(3-17b) 

and  the  expansion  of  g(cP.<P,)  is  performed 

as  follows.  Using  the  relations^” 

me  («;h2)  -fl  ce  (z;h2) 
m  m 

m  =  0,  1,  2  •  •  • 

(3-  18a) 

mejz;h2)  =i JZ  se  (z;h2) 

-  m  m 

m  =  1 , 2, 3  ... 

(3-18b) 

(z;h)  -  Mc^  (z;h) 
m  m 

m  =  0,  1, 2  ... 

(3-18c) 

( z;h)  =  (-l)mK&J(z;h) 

-m  m 

m  =  1, 2, 3  ... 

(3-18d) 

one  may  write 

00  00 

g((P.cp,)  =  2i  gmte(P')  =-4i  ^  { - l)m  se 
m=l  m=l 

(n')se  (-t))  - rxr - 

m  m  Ms  3  (0) 

m 

(3-19) 

because 

Mc(1)  (0)  =  0  for  m  =  0,  1,2  •*  . 
m 

(3-19a) 

Furthermore 

^  fuolh2)] 

Ms\3,{0)  2  L  j 

(3-20a) 

Mi<2>(0>  -in .4  f  ^  2,1 

^7n^-TTh  Hh)J 

(3 -20b) 
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se^Tijh2)  =  sin  cp  -  sin  +  0{h4)  (3-20c) 

2  h2  4 

sc^rp  h  )  =  sin2cp--~  sin4Cp  +  0{h  )  .  (3~20d) 

With  Eqs.  (3- 19),  ( 3-20)  and  similar  expressions  for  m>2  it  can  be  shown  that 
gm(<p.cp')  «  0(^)m 
The  leading  term  has  the  form 

g^cp.cp')  =  4nh2  sincp  sincp'  +  0(h4)  .  (3-21) 

2 

The  quantities  cp,<p'  and  h  may  be  expressed  in  terms  of  0,9',  a,  etc. 

Rather  than  using  the  general  Eqs.  (2-6),  (3-3),  (3-5),  (3-7)  and  (3-10a),  one 
may  obtain  the  result  by  inspection  for  the  special  case  under  discussion  (see 
Fig.  (5)  ): 


Fig.  5  The  conducting  strip  in  the  real  (x,  y) 
and  the  transformed  (u,  v)  space. 


22 


R 

■  sin  6 

cos  6 

sin  <p 


^  =Ju2  +  V2  = 

JZ.  4  Z  =  rN<6> 

^  5  Ja? 

(3-22a) 

< 

o 

1 

o 

1 

JIT  cos  8q 

(3-22b) 

Ju2  +  v*  r  N(8  ) 

1  o  o  o  o' 

N(8o) 

"o  .  -f*  O 

Jisin  0q 

(3-22c) 

Ju^  +  v^  r  N(8  ) 

V  o  o  o  o 

N(0q) 

V  _  -usinfi  +vcos6 

-xcoe  0  +  y  sin  8 

O  1  O  r - 

7WI  N{9o) 


tae  = 


-sin0  cos  8  +  cosS  sin0  sin (6  -6) 

- N(0  ) - -£7^  -Fpr°^Ty 

o  o 


(3-22d) 


Thus,  Eqs.  (3- 17a),  (3-21)  and  (3-22d)  yield 

7  sin(0  -0')  sin(8  -8)  . 

fil^co.cp')  sgj(9f  0')  =  TT(kQa)  ae  - —  +  0(koa)4  (3-23) 


One  may  verify  that  Eq.  (3-23)  represents  the  radiation  from  an  electric  dipole 
line  source  and  find  therefrom  the  polarizability  of  the  strip. 


D.  ELECTRIC  DIPOLE  LINE  RADIATION,  AND  THE  POLARIZABILITY 
OF  A  STRIP 

Consider  Maxwell's  equations 

7xE  =  ium^  (3 -24a) 

7xH  =  iir  e  •  E  +  J  (3-24b) 


with  e 
in  the 


given  by  E 
z-direction 


(2- la),  H  given  by  Eq,  (2-3b)  because  there  is  no  variation 
">r —  «  o]  and  the  source  term  given  by 


I 
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5  ( r) 


£  =  (j^«in0o  +  ioco8eo)  p 


(3-24c) 

(3-24d) 


This  excitation  describes  a  uniform  line  distribution  of  electric  dipoles  oriented 
perpendicular  to  the  line  axis.  From  Eqs.  (3-24),  one  may  derive  the  wave 
equation 


[* 


a2  a2 

— T  +  e  — 7  + 
Sx*  dy^ 


3J  9J 

k  J  H  =  e  ' a  =  A6  &  » 


where  A  =  p 


e  sin  a 


-a  cos  0 


o  -s] 


o  57 

If  we  define  H=-  AG,  then  the  equation  for  G  reduces  to 

a  e^~T  +  k2G  =-  6(r) 

?x  dy 


(3-25) 


(3-25a) 


(3-26) 


whose  solution,  subject  to  a  radiation  condition  at  infinity,  has  already  been 
given  by  Gd(jr,  0)  in  Eq.  (3-14).  Thus,  the  solution  of  Eq.  (3-25)  is 


H(r)=lp[ 

=  ipkoae 


a  cos  8_  15—  -s  sin  0Q  -gp  j  H'^ 


d 

o  ^ 


cos  9 


sin  9 


.  sin(9-9  ) 

=TP  ko  a€ 


^ay2+  ex  2  ^ 

NT&r  +  ex2j 


ay  +  cx 

I— 


2  2 
ay  +cx 


For  kJay*+  cx2  =  k  rN(9)»  1  one  may  write 


[korN(0)  -  -i  H(^  JkorN(9)J 
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so  that 


H  (r)  ~  pk  ae 


sin  (0  -  eo)  t  . 


N  (0) 


H' 


4  o 


(1) 


korN(0) 


l 


comparison  of  Eqs.(3-15),  (3-27)  and  (3-23)  yields 

.2 


P  = 


TT 


(koa)‘ 


sin(0  -07) 

V  t  A\/l-  _  \  * 


— 2__  +  0(k  a) 
N  (0  ')  o 


(3-27) 


(3-28) 


where  Aq  is  the  intensity  of  the  incident  field  at  the  location  of  the  obstacle. 
Eq.  (3-27)  shows  the  leading  term  in  the  expansion  of  the  scattered  field  as  an 
electric  dipole  term,  and  Eq.  (3-28)  gives  the  intensity  of  the  dipole  induced 
in  the  strip  by  the  incident  field.  The  first  order  approximation  of  the  scat¬ 
tering  pattern  as  given  by  Eq.  (3-23)  is  shown  in  Fig.  6.  A  calculation  simi¬ 
lar  to  the  above  for  the  determination  of  the  quadrupole  radiation  is  given  in 
appendix  B. 


Fig.  6  Dipole  radiation:  first  order  approximations  of 

the  scattered  field. 
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E.  CONCLUDING  REMARKS  ON  LONG-WAVE  SCATTERING 

The  expansion  of  Eq.  (3-19)  can  be  carried  out  to  higher  power*  of 

ka  (or  h).  In  the  same  manner  as  in  the  preceding  section,  one  may  show 

4 

that  the  leading  term  in  g^cp.cp  )  is  0(kQa)  and  its  angular' pattern  is  that 
of  a  linear  electric  quadrupole. 

If  one  considers  the  problem  of  an  elliptic  cylindrical  scatterer  (b£0), 
the  leading  term  in  the  exapnsion  of  the  scattered  field  is  from  the  m  =  0 
term  in  Eq.  (3-15a).  This  respresents  an  induced  "monopole*,  i.  e. ,  an 
equivalent  line  source.  There  will  also  be  two  Induced  dipole  terms,  one  along 
the  major  axis  and  one  along  the  minor  axis  of  the  elliptic  scatterer  ,  and  Eq. 
(3-28)  will  e  written  in  dyadic  form 

p  s  P  •  E 

(12) 

where  E  is  the  incident  electric  field.  The  expansion  has  been  performed  1 
up  to  O(h^).  By  properly  transforming  sin  mep  and  cos  m  cp  (as  has  been  done 
with  sin  cp  in  Eqs.  (3-22)),  one  may  find  the  scattering  pattern  in  the  anistropic 
medium. 

For  the  case  where  the  uniaxial  medium  represents  a  plasma  in  a  strong 
constant  magnetic  field,  a  is  real  and  positive,  but  e  may  be  positive  or 
negative  (in  the  lossless  case)  or  complex  with  positive  imaginary  and  positive 
or  negative  real  part  (in  the  lossy  case).  In  the  former  instance,  it  is  seen 
from  Eq.  (2-11)  that  if  e  <  O,  N(0)  will  be  imaginary  in  that  part  of  space  for 
which 


tan  G  >  tan  9 


(3-2  <* 


Felsen  ^  has  shown  that  if  Im  N(8)  >  O  (for  the  e-  time  dependence),  the 
expression  given  by  Eq.  (3-13  )  for  the  field  radiated  by  a  line  source  in  an 
infinite  medium  (which  is  identical  with  G^  of  Eq.  (3-11)  )  still  satisfies  the 
radiation  condition.  We  may  define  the  square  root  in  Eq.  (2-11)  to  be  positive 
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when  real  and  have  a  positive  imaginary  part  when  complex.  This  restricts  the 
s 

quotient  —  to  have  a  positive  imaginary  part, 

lm(|)>0  (3-30) 

It  i?  *■<">  «how  that  the  same  condition  on  —  is  necessary  for  the 

scattered  part  G  of  Eq.  (3-11)  to  be  convergent  and  to  satisfy  the  radiation 

condition.  From  theorem  1  in  section  1.  3  of  reference  (8)  one  may  deduce  that 

the  various  Mathieu  functions  and  their  derivatives,  which  appear  in  Eq.  (3-11), 

are  analytic  functions  of  h  for  fixed  (possibly  complex)  arguments.  Thus,  for 

fixed  r,  r'  and  ?  ,  G  is  an  analytic  function  of  h.  It  follows,  that  if  the 
-  os 

power  series  expansion  of  G  converges  within  a  certain  real  interval 
-hQ<  h<  hQ,  it  will  converge  within  the  circle  |h|<|hoi  for  complex  values 
h*.  Also  (according  to  reference  (8)  p.  98)  the  various  formul:-  s  for  the  Mathieu 
functions  of  which  we  make  use  in  this  chapter,  are  valid  for  arbitrary  complex 
arguments  and  parameters.  Thus  Eq.  (3-15)  is  valid  even  if  R,  R/,cp  and  cp7 
are  complex.  However,  in  order  to  satisfy  the  radiation  condition,  the  inequality 
(3-30)  has  to  be  satisfied,  because  in  (3-15)  G  is  seen  to  be  asymptotically 
proportional  to  G^  (multiplied  by  the  scattering  pattern).  The  same  singulari¬ 
ties  which  appear  along  the  direction  0,  (Eq.  (3-29)  )  in  the  direct  field  will 

C 

appear  in  the  scattered  field  for  negative  real  —  .  Eq.  (3-28)  also  shows  that  the 
induced  dipole  (and  higher  multipoles)  would  be  inflate  if  8/  =  0  ,  because  of  the 
singularity  of  the  incident  field  along  that  direction. 

The  scattering  pattern  function  g  (Eq.  (3-15)  )  in  isotropic  media  is  a 
function  of  the  geometry  of  the  scatterer,  and  its  orientation  with  respect  to  the 
directions  of  incidence  6',  and  observation  9.  Eq.  (3-23)  shows,  that  the 
orientation  with  respect  to  the  optic  axis  9q  is  an  additional  parameter  in  aniso¬ 
tropic  media.  Also  the  ray  refractive  index  in  the  directions  0  and  0*  enters 

explicitly  into  the  expansion  of  g.  The  ray  refractive  index  is  a  property  of  the 

(13) 

anisotropic  medium  ,  and  may  be  found  in  more  general  anisotropic  media  as 
well.  Therefore  the  results  obtained  indicate  what  properties  of  the  medium  enter 
the  expression  of  the  scattering  pattern  g  or  the  scattering  cross-section  which 
is  related  thereto  when  the  scatterer  is  embedded  in  an  anisotropic  medium. 

*See  for  example  chapter  9  in  "Advanced  Calculus"  by  W.  Kaplan  (Addison- 
Wesley,  1952)  for  proofs  of  these  statements. 
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CHAPTER  IV 

DIFFRACTION  IN  THE  SHORT  WAVE  RANGE 


A.  INTRODUCTION 

Diffraction  by  objects  of  relatively  arbitrary  shape  is  known  to  arise 

from  certain  localized  regions  on  the  scatterer  surface  when  the  obstacle 

dimensions  are  large  compared  to  the  wave  length  of  the  incident  radiation. 

The  resulting  phenomena  are  described  conveniently  in  terms  of  rays  which 

account  for  diffraction  effects  in  addition  to  those  of  reflection  and  refraction 

in  geometrical  optics.  For  electromagnetic  wave  propagation  in  isotropic 
(3) 

regions,  Keller'  has  presented  a  ray  theory  for  the  construction  of  the  geo¬ 
metric-optical  as  well  as  the  diffracted  field.  The  theory  is  based  on  certain 
postulates  which  have  been  verified  in  special  but  representative  caoes  by  com¬ 
parison  with  the  asymptotic  representation  of  rigorous  solutions.  The  purpose 
of  the  present  investigation  is  to  furnish  a  similar  interpretation  when  the 
medium  surrounding  the  obstacle  is  anisotropic. 

A  significant  difference  between  propagation  in  isotropic  and  anisotropic 
regions  is  the  distinction  in  the  latter  between  the  directions  of  propagation  of 
the  phase  fronts  (wave  normal  direction)  and  of  the  energy  (ray  direction)  in  a 
plane  wave.  Since  the  fields  along  a  ray  are  locally  those  of  a  plane  wave,  die 
laws  of  reflection  of  a  ray  may  be  deduced  from  a  solution  of  the  corresponding 
plane  wave  boundary  value  problem.  Alternatively,  these  laws  may  be  deriveu 
by  applying  Fermat’s  principle  in  a  form  suited  to  propagation  in  anisotropic 
regions,  and  this  approach  is  adopted  here.  The  field  amplitudes  along  a 
reflected  ray  may  be  determined  from  the  principle  of  conservation  of  energy 
in  a  ray  tube,  thereby  permitting  the  construction  of  the  reflected  field  by 
arguments  of  geometrical  optics.  Diffraction  effects  may  arise  from  surface 
singularities  or  from  the  vicinity  of  the  shadow  boundary  on  the  obstacle.  Two 
types  of  representative  problems  are  investigated  in  this  connection:  diffraction 
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by  a  perfectly  conducting  half  plane,  and  diffraction  by  smoothly  curved  sur¬ 
faces.  In  each  case,  the  solution  is  obtained  from  the  asymptotic  evaluatior  of 
a  rigorous  formula  and  its  subsequent  interpretation  in  general  ray-optical 

terms,  thereby  permitting  a  generalization  to  anisotropic  regions  of  the  geo- 

(3) 

metrical  theory  of  diffraction  .  The  formulation  exhibits  explicitly  the  depen¬ 
dence  on  the  local  properties  of  the  medium  and  the  scatterer  surface. 

A#  meaiioneu  in  Chapter  II,  the  present  investigation  is  restricted  to 
perfectly  conducting  cylindrical  structures  excited  by  axially  independent  inci¬ 
dent  fields,  and  most  of  the  phenomena  described  pertain  to  uniaxial  anisotropy, 
with  the  optic  axis  of  the  medium  perpendicular  to  that  of  the  scatterer.  The 

diffraction  problem  may  be  reduced  by  Felsen's  method  of  scaling  to  an  equiva- 

(2) 

lent  one  in  an  isotropic  region'  ,  thereby  permitting  the  direct  construction  of 
rigorous  asymptotic  solutions  and  their  ray-optical  interpretation.  The  analytic 
continuation  of  the  solution  from  arg  e  =  0  to  0  <  arg  e  <  tt  is  given  specific 
attention.  The  continuation  is  justified  in  detail  for  the  half  plane  problem  and 
the  parabolic  cylinder  problem,  thereby  lending  an  extended  range  of  validity 
to  the  results  derived  by  the  scaling  technique. 

B.  GEOMETRICAL  OPTICS 

It  is  well  known  that  geometrical  optics  predicts  correctly  the  dominant 
effects  of  the  electromagnetic  field  when  the  wave  length  tends  to  zero.  The 
geometric  optical  field  is  comprised  of  the  incident  (direct),  reflected  and 
refracted  constituents  and  its  properties  are  described  conveniently  in  terms  of 
rays.  Since  the  geometric  optical  field  is  locally  a  plane  wave  field,  the 
reflection  and  refraction  properties  of  the  rays  are  deducible  from  those  of 
the  corresponding  plane  waves.  Moreover,  the  local  reflection  and  refraction 
characteristics  of  slowly  curved  interfaces  are,  in  the  geometric-optical 
approximation,  the  same  as  for  a  plane  interface  tangent  at  the  point  of  i  npact 
of  the  ray,  so  that  the  information  extracted  from  the  analysis  of  plane  wave 
reflection  and  refraction  by  a  plane  interface  separating  two  media  suffices  for 
a  determination  of  the  initial  direction  and  amplitude  along  a  reflected  or 
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refracted  ray.  It  must  be  kept  in  mind  that  in  an  anisotropic  medium,  the  ray 

direction  (direction  of  local  energy  transport)  differs  from  that  of  the  wave 

normal  (direction  of  progressing  phase  fronts),  with  the  relative  orientation  of 

the  ray  and  wave  normal  vectors  determined  by  the  medium  constants.  The 

utility  of  the  refractive  index  surfaces  for  the  anisotropic  medium  in  predicting 

the  plane  wave  (and  therefore  the  ray)  propagation,  reflection  and  refraction 

I  K\/A\ 

characteristics  has  been  empha.sized  elsewhere'-'*'''"'.  In  the  present  analysis, 
the  ray  trajectories  are  derived  from  purely  ray-optical  considerations  via 
Fermat's  principle  of  stationary  propagation  time,  and  it  is  shown  that  the 
results  so  obtained  agree  with  those  derived  from  the  previously  mentioned 
plane  wave  considerations.  Since  only  perfectly  conducting  surfaces  are  con¬ 
sidered,  the  analysis  involves  the  direct  and  reflected  rays  only. 

The  pertinent  form  of  Fermat's  principle  for  propagation  along  a  ray 
path  s  in  an  anisotropic  medium  is  (reference  6,  p.  289), 

ftj  Ndl  =  0  (4-1) 

s 

where  N,  the  ray  refractive  index,  is  related  to  the  ordinary  refractive  index 
n  via  N  =  n  cos  y»  with  Y  denoting  the  angle  between  the  ray  and  the  wave 
normal.  Since  the  ray  refractive  index  determines  the  propagation  speed  of  the 
phase  front  along  the  direction  of  the  ray,  the  "stationary  time"  in  Fermat's 
principle  refers  to  the  latter  and  not  to  the  time  of  energy  transport.  The 
dependence  of  N  on  n  may  also  be  expressed  as  (reference  6,  p.  253) 


N 

n  cos  Y  , 

(4- 2a) 

tan  Y  = 

1  d  n  (cp) 

n(cp)  &p  * 

(4- 2b) 

cp 

a  -  y  » 

(4-2c) 

where  the  presence  of  a  variation  of  n  with  direction  is  indicative  of  the 
medium  anisotropy.  The  singles  9  and  cp  measured  from  the  y-axis  identify 
the  directions  of  the  ray  and  wave  normal,  respectively.  (Fig.  7) 
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For  the  class  of  problems  considered  here  the  refractive  index  may  be 
calculated  by  assuming  a  j  wave  solution  of  the  form 


A 

,  ik  nk .  r 

H(x,  y)  =  Ae1  x  +  *  *>  =  Ae't-i*  Ae  0  ~  ", 


(4-3) 


substituting  into  the  homogeneous  Eq.  (2-4),  and  determining  thereby  the  plane 
wave  dispersion  relation  which  connects  the  wave  numbers  £  and  n  : 

K  =  -  (kQ  A  -  aC  )  (4-4) 

with  A  defined  by  Eq.  (2-4b).  In  Eq.  (4-3), 

k  =  x  C+V  h  =  k  nfc  (4-5) 

—  — o  o  o  — 


is  the  wave  normal  vector  defining  the  direction  of  propagation  of  the  phase 

A 

fronts,  k  =  k/k  is  a  unit  vector  in  the  direction  of  k,  and  _r  =  xox+^oy 
is  the  position  vector.  As  noted  previously,  the  energy  flux  vector  S  (which, 
in  a  lossless  medium,  is  the  time-averaged  Poynting  vector  Re  JE  x  H*)  is 
inclined  with  respect  to  k  by  the  angle  y,  and  the  various  angles  employed 
subsequently  are  schematized  in  Fig.  7.  From  Eqs.(4-4)  and  (4-5), 


n(cp) 


l 


2  .  2 
€  cos  cp  +  a  sin  cp 


(4-6) 


and  the  above-mentioned  refractive  index  diagram  is  obtained  by  plotting 
n(cp)  vs.  cp.  If  the  ray  progresses  along  the  9-  direction,  with  tan0  =  x/y, 

A 

then  (£x  +  ny)  =  kQn(cp)k.£  =  kQrN(9),  and  it  is  not  difficult  to  show  from 
these  considerations  or  from  Eqs.  (4-2)  that 


(4-7) 


which  reduces  to  the  expression  given  by  Eq.  (2-11)  when  8=0  in  Eq.  (2- lb). 


Fig.  7:  Plane  wave  propagation  in  an  anisotropic  medium. 


It  has  been  mentioned  in  Chapter  II  that  two  special  cases  arise  which 
are  of  interest  in  connection  with  plasma  media  subjected  to  an  external  steady 
magnetic  field.  If  the  magnetic  field  is  parallel  to  the  z-axis,  then  a  = 
so  that  both  n(c£)  and  N(3)  are  constant  and  equal  to  A.  The  wave  equation 
(2-4)  then  reduces  to  the  one  for  an  isotropic  medium,  and  the  anisotropy 
manifests  itself  only  through  ‘he  relations  (2-5  a,  b)  which  differ  from  the 
isotropic  case  in  view  of  the  non-vanishing  3.  Since  the  anisotropic  features 
are  caused  in  this  instance  primarily  by  the  presence  of  boundaries  rather 


than  by  the  medium  itself,  this  case  is  not  considered  further.  Instead,  we 

2 

specialize  to  the  uniaxial  situation  3  =0,  a  =  1,  e  =  1  -  (Wp/io)  ,  which 
arises  when  the  external  magnetic  field  is  very  strong  andi  is  oriented  along 


one  of  the  coordinate  axes  perpendicular  to  z,  say  y  (uu^  denotes  the  plasma 
frequency  based  on  the  motion  of  electrons  only,  and  u>  is  the  frequency  of 


the  electromagnetic  field).  Eq.  (4-7)  then  reduces  to 
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N{6)  =  Vcos^0+esin8  ,  8  =  0,  a  =  1 


(4-8) 


With  the  ray  refractive  index  specified,  the  ray  trajectories  may  be 
calculated  from  Eq.  (4-1).  Since  the  medium  is  assumed  to  be  homogeneous, 
the  rays  proceed  along  straight  lines  and  it  is  trivial  tc  determine  the  direct 
path  from  the  source  to  the  observation  point.  If  a  plane  reflecting  surface  is 
present  as  in  Fig.  8  inclined  at  an  angle  0Q  to  the  "optic  axis"  y,  then  the 
optical  path  from  point  P^.yj)  to  point  Qfx^,  y2), 

Q  _ _ 

L  =  J  Ndl  =l/(y1-y)^  +  efxj-x)2  +  4  (y-y2)2  +  e(x-x2)2  (4-9) 


must  be  extremized  subject  to  the  constraint  that  point  C  on  the  trajectory 
lies  on  the  plane. 


Fig.  8:  Reflection  of  a  ray  from  a  perfectly  conducting  plane 
in  an  anisotropic  medium. 
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The  details  of  the  calculation  are  given  in  Appendix  C  and  the  result  furnishes 
a  relation  between  the  reflected  ray  direction  8  ,  the  incident  ray  direction 
9.,  and  the  angle  0  ,  all  angles  being  measured  in  the  clockwise  sense  from 
the  y-axis: 

-  ( 1  -  e  tan^  0  )  tan  0.  +  2  tan  0 

tan  0  =  - * - - - i - —  (4-10) 

(1-etan  8  )+2etan8  tan 8. 

o  o  i 


This  foi-mula  may  be  shown  to  agree  with  Eq.  (10)  of  reference  15  which  was 
obtained  from  a  solution  of  the  boundary  value  problem.  If  all  angles  are 
measured  from  the  normal  to  the  plane,  the  formula  becomes  instead 


2A  (8  ) 

ta„ar  = 


(4- 10a) 


where  A_(0 I  )  and  A?(9  )  are  defined  by  Eqs.  (4-24),  8  being  the  angle 

between  the  conducting  surface  and  the  optic  axis,  measured  clockwise  from 
this  axis.  In  the  isotropic  limit  e  =  1,  this  expression  yields  correctly  the 
specular  condition  tan  =  tan  a..  The  preceding  results  may  also  be  applied 

to  a  smoothly  curved  surface  since  reflection  takes  place  locally  as  on  an 
infinite  tangent  plane.  With  the  ray  trajectories  known,  it  follows  from  the 
definition  of  the  ray  refractive  index  that  the  phase  change  over  a  distance  d 
along  a  ray  directed  at  an  angle  9  with  the  y-axis  is  kQdN(0).  Referring  to 
Fig.  8,  for  example,  the  phase  increment  along  the  incident  ray  over  the  dis¬ 
tance  PC  between  points  P  and  C  is  kQ  PC  N(8j),  while  t>  j  corresponding 
quantity  along  the  reflected  ray  between  points  C  and  Q  is  k  CO  N(8?). 

Attention  may  now  be  given  to  the  field  amplitudes  along  the  direct  and 
reflected  rays.  If  the  excitation  is  in  the  form  of  a  line  source  of  magnetic 
current  located  at  the  point  (x/,y/),  the  magnetic  field  H(x,  y)  is  given  by  the 
solution  of  Eq.  (2-4)  subject  to  a  radiation  condition  at  infinity.  The  solution  may 
be  constructed  by  Felse,n's  method  which  reduces  the  differential  operator  :in 
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Eq.  (2-4)  to  the  ordinary  Laplacian  in  the  u-v  space.  The  corresponding 
Green's  function  is  known  and  may  then  be  transformed  back  into  the  x-y  frame, 
with  the  result  (reference  2,  Eq.  (28): 


H(x,  y) 


a1  e  A 
o 


H  <» 


(4-11) 


which  agrees  with  Eq.  (3-13).  The  energy  flux  S  =  Re  (IS  x  H*)  is  calculated 


from  Eqs.  (2-5)  and  (4-11)  and  yields  the  asymptotic  result  (when 
k  N(0)  |r  -  r#|  »  1), 


S 


r 

— o 


TW  1  H<x»  y)l 


(4-12) 


thereby  confirming  the  straight-line  character  of  the  rays,  and  showing,  in 
addition,  that  the  rays  emanate  radially  from  the  source. 

While  Eqs.  (4-11)  and  (4-12)  have  been  derived  here  on  the  assumption 

that  a  and  e  are  positive,  it  has  been  shown^  (for  the  uniaxial  case  8=0,  a=l,) 

that  Eq.  (4-11)  remains  valid  also  for  0<arge<n  provided  that  the  resulting 

N(6)  is  defined  to  have  a  positive  imaginary  part.  The  latter  requirement  is  in 

accord  with  the  radiation  condition  and  assures  the  decay  of  the  fields  at  infinity, 

either  due  to  the  presence  of  dissipation  when  €  is  complex  or  due  to  shadow 

effects  when  s  is  negative  real.  When  s<0,  (for  the  uniaxial  plasma,  this 

occurs  when  uKju  )  S  in  Eq.  (4-12)  differs  from  zero  only  in  those  angular 
P  2  21/2 

regions  wherein  N(6)  =  (cos  9  +  ssin  S)  is  real.  Illumination  is  there¬ 
fore  confined  to  a  wedge-shaped  region  centered  at  the  source  and  surrounding 
the  optic  axis.  A  simple  ray-optical  description  obtains  when  the  Hankel  func¬ 
tion  in  Eq.  (4-11)  can  be  replaced  by  its  asymptotic  form,  i.  e. ,  when 
k  N(8)  |r  -  r'|  »  1.  While  this  condition  can  always  be  realized  for  sufficiently 
large  kQ  when  N(0)  ^  0,  it  fails  on  the  shadow  boundary  (N(@c)  =  o).  In  the 
transition  region  surrounding  @c  =  tan”*  (1/^Ul ),  one  must  employ  the  exact 
formula  (4- 11). 
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Away  from  the  transition  region  svir rounding  0^,  the  asymptotic  form  of 
Eq.  (4-li)  maybe  employed  to  furnish  the  magnitude  of  H  at  a  distance  Rj 
from  the  source,  and  the  corresponding  flux  density  at  (x^.y^)  is  then 
obtained  from  Eq.  (4-12).  One  may  now  apply  the  principle  of  conservation  of 
energy  to  a  narrow  tube  of  rays  to  deduce  that  dA^  =  f>2  dA^,  where  dA^  i* 
the  cross- sectional  area  of  the  iube  c*  Rj,  and  dA2  i*  i)»e  analogous  quantity  at 
a  distance  R2>Rj  along  the  ray.  For  this  two-dimensional  configuration, 
dAj  =  Rj  d6,  dA^  =  R^  d9,  so  that 

|H(x2,  y2)|  =  iHlx^yj)!.^^  .  (4-13) 

When  an  incident  ray  strikes  a  smoothly  curved  surface  as  in  Fig.  9,  it 
gives  rise  to  a  reflected  ray  whose  direction  may  be  inferred  from  the  reflection 
law  in  Eq.  (4-10),  with  9q  representing  the  angle  between  the  optic  (y)  axis 
and  the  tangent  to  the  surface  at  the  point  of  reflection. 


Fig.  9:  Geometric -optical  construction  of  the  field 
due  to  a  line  source  in  the  presence  of  a  smoothly  curved 
perfectly  conducting  cylinder. 
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To  determine  the  reflected  ray  amplitude,  the  conservation  of  energy  argument 
is  applied  to  a  narrow  tube  of  reflected  rays  which  appears  to  emanate  from  a 
focus  F  located  at  a  distance  "a"  behind  the  reflecting  surface.  Upon 
applying  the  sine  law  to  the  three  triangles  SAC,  FAC  and  RAC  in  Fig.  9 
one  obtains  the  relations 


ds 

- 

o 


b 


(4-1 4a) 


ds  d' 

d0 1  sin(8j-So) 


(4-14b) 


ds  _  ds  _  a 

d$-,  m  d6 .  +  n  d9  sin  (8  -6,) 

2  1  o  o  2 


(4-14  c) 


where  we  have  utilized 


d9,  =  m  d8 ,  +  n  d8  , 

2  1  o 


(4-14d) 


and  the  partial  derivatives  m=  and  n^dO^/dO^)  can  be  evaluated 

from  Eq.  (4-10);  b  is  the  radius  of  curvature  of  the  cylinder  at  the  reflection 
point;  "a"  can  now  be  found  as  follows: 


l 

a 


m 

d7 


sin(9, -0  ) 
1  o 

sin<V02J 


sin(5o-  62) 


(4-15) 


For  the  special  case  of  an  isotropic  medium,  6^  =  20^-0  j  and  Eq.  (4-15) 
reduces  tc  the  known  formula  (of  reference  16,  Eq.  (3)) 


a 


1 

d7 


+ 


2 

b  cos  a 


(4- 15a) 


(see  F’ig.  9) 
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The  required  geometrical  quantities  are  now  determined  and  we  may 
proceed  to  the  consideration  of  the  power  flow  in  the  incident  and  reflected  ray 
tubes.  If 


H 

n 


-«o  v*  / - 2 -  -in/4 

- 4 —  J  ~~Tjnrr  e 

»  "o  '  1' 


(4-16) 


denotes  the  incident  field  at  unit  distance  from  the  source  (see  Eq.  (4-11),  then 

the  incident  field  H.#  at  the  reflection  point  B  in  Fig.  9  is  given  via  Eq.  (4-13) 

i  - 1  /■* 

and  the  inclusion  of  the  phase  change  as  H.'  =  H  (d7)  ”  **  exp  (ikQ  d7  N(9j) ); 

the  associated  energy  flux  in  the  incident  ray  tube  is 

S.  dA.  =  7i  i  H  | 2  fd'NO.)]'1  d7d0. 
ii  o  o’  i  l'j  l 

(see  Eq.  (4-12)).  Similarly,  if  denotes  the  reflected  magnetic  field  at 
B,  then  the  flux  in  the  reflected  ray  tube  at  the  surface  is 

SrdAr  *  rolHrl2[N(52>]'1  adS2  • 

By  conservation  of  power,  S^dA^  =  S.  dA. ,  thereby  permitting  the  determina¬ 
tion  of  |H^|  in  terms  of  i  H.7 1  and  the  derivative  (dS^/dS.^)  in  Eqs.  (4-14  b.  c). 
|Hrl  at  the  observation  point  P  in  Fig.  8  is  then  given  by 


and  the  phase  of  differs  from  that  of  by  kodN(8^).  Since  the  plane 
wave  reflection  coefficient  for  the  magnetic  field  is  easily  shown  to  be  equal  to 
(+1),  no  phase  change  occurs  upon  reflection,  and  one  finds  by  combining  the 
various  expressions  that  the  reflected  field  at  P  due  to  excitation  by  a 

line  source  at  S  is  given  by: 


H  = 
r 


J  tb‘T-  / 

2V2rT  'y]  N^Qj)  sin(0o-8. 


i^JdjNtGj)  +  d  N(82)J-i  tt/4 


dmsintB^  80)"  h'  d“ 

sin(0Q  -  to?)  b  sin(0 


(4-17) 


The  incident  field  IT  over  the  direct  path  D  from  S  to  P  is  found  from  Eqs. 
'4-16)  and  (4- 13)  as: 


H.  = 


-'JU  £  Jz 

o4 


1  2^2  n  N(83)koD 


ik  DN(0_)  -  in/4 
e  o  5 


(4-18) 


The  geometric-optical  field  H  =  H. +  Hy  may  be  expected  to  furnish  good 
approximation  to  the  actual  high-frequency  field  in  the  illuminated  region  of  the 
smoothly  curved  scatterer.  This  statement  is  verified  subsequently  from  a 
study  of  the  asymptotic  field  solution  for  the  special  case  of  a  parabolic  cylinder. 

It  may  also  be  shown  that  the  geometrical  optics  formula  in  Eq.  (4-17) 

yields  the  correct  result  when  the  medium  is  taken  as  isotropic  (n(0)  -  l),  or 

when  the  reflecting  surface  is  planar  (b=t1).  In  the  former  case,  the  relation 

between  9^,  8j  and  0Q  simplifies  to  02  =  20  -9  Upon  introducing  tv.  angle 

a  =  (tt/2)  -  (0,-8  )  between  the  incident  ray  and  the  surface  normal  (see  Fig. 

1  o  (171 

9)  one  may  reduce  Eq.  (4-17)  to  the  known  formula  '  ”  ' 


i  k  (d*  +  d)  -  i  tt/4 
o 


r  zjrn  jn 

d  +  d'  +  — 

*v  0 1 

b  cos  a 

(4-19) 


Alternatively,  when  b  =  ",  Eq.  (4-17)  agrees  withthe  expression  for  reflection 
of  a  cylindrical  wave  from  an  infinite  plane,  which  was  obtained  from  a  solution 
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(15) 

of  the  wave  equation  .  This  statement  is  verified  in  Appendix  D.  Since  the 
latter  result  was  shown  to  apply  as  well  when  0<  arg  e  <  n,  the  same  range  of 
applicability  may  be  expected  from  Eq.  (4-17)  provided  that  the  imaginary  part 
of  N(0)  is  restricted  to  be  positive.  If  e  <  0,  the  radiation  emanating  from  the 
line  source  is  confined  to  a  wedge  shaped  region  ]  tan  6j  |<  |  tan  9  |  and  the 
simple  ray-optical  description  fails  at  the  shadow  boundaries  |  tan  0^  |  =  tan  0c> 
thereby  necessitating  use  of  the  exact  formula  (4-11).  Analogous  considerations 
are  val  r  for  the  reflected  field  which  is  also  confined  to  the  angular  interval 
|  tan  0~|  <  I  tan  0  J  ;  since  the  reflected  field  appears  to  come  from  an  image 

w  C 

source,  it  is  to  be  expected  that  its  value  near  the  shadow  boundary 
|  tan  0^1  =  tan  may  be  obtained  by  retention  of  the  Hankel  function  in  (4-11) 
with  the  appropriate  argument. 

C.  DIFFRACTION  BY  A  STRAIGHT  EDGE 

The  geometric -optical  solution  constructed  in  the  preceding  section  may 
be  expected  to  yield  the  dominant  contribution  in  the  illuminated  region  of  the 
obstacle,  but  it  yields  no  information  about  diffraction  effects  which  modify  the 
field  in  this  region  and  ,'J.so  account  for  field  penetration  into  the  shadow  zone. 
Diffraction  phenomena  arise  when  the  surface  contains  singularities  such  as 
edges  or  corners,  and  also  from  the  vicinity  of  the  shadow  boundary  on  a 
smoothly  curved  object.  A  prototype  structure  lor  the  study  of  edge  effects,  a 
perfectly  conducting  half  plane,  is  investigated  in  this  Section.  When  the  medium 
is  uniaxial  (6  =  0,  a  =  1  in  Eq.  (2-lb)),  the  previously  mentioned  method  of 
scaling  may  be  employed  to  reduce  the  problem  of  scattering  by  perfectly  con¬ 
ducting  cylindrical  surfaces  in  the  anisotropic  medium  to  equivalent  conventional 

(2) 

problems  in  an  isotropic  region.  The  latter  configuration  is  a  real  one  when 
e  is  positive  real,  but  loses  its  physical  significance  for  other  values  of  €.  In 
utilizing  the  auxiliary  isotropic  space,  it  is  therefore  necessary  to  construct  the 
solution  for  e>  0  ar.J  to  seek  an  extended  range  of  validity  by  analytic  con¬ 
tinuation. 

The  problem  of  diffraction  by  a  perfectly  conducting  half  plane  in  a 
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(2) 

uniaxially  anisotropic  region  has  already  been  treated  from  this  viewpoint 

and  the  solution  has  been  given  in  a  form  analogous  to  that  derived  first  by 

(18) 

MacDonald. '  '  In  this  representation,  it  is  convenient  to  employ  angles 

which  delimit  the  various  ray-optical  domains;  while  the  interpretation  is 
unambiguous  when  e  >  0,  it  is  not  clear  how  to  effect  the  continuation  into  the 
range  of  complex  or  negative  real  e  since  the  angles  then  become  complex. 

To  clarify  the  picture,  an  alternative  analysis  has  been  carried  out 
wherein  distances  rather  than  angles  are  employed  throughout.  The  problem 
has  been  solved  directly  by  the  Wiener-Hopf  technique  without  the  intervention 
of  the  scaling  argument,  and  the  solution  cast  into  a  form  which  facilitates  the 
analytic  continuation.  The  same  basic  representation  may  also  be  obtained  when 

Felsen's  method  is  applied  judiciously  to  the  Wiener-Hopf  solution  in  an  iso- 

(4)  (19) 

tropic  medium  .  This  approach  has  been  adopted  elsewhere  ,  and  serves  as 

a  check  on  the  correctness  of  the  direct  approach. 

The  physical  configuration  is  shown  in  Fig.  10,  with  the  half  plane 

occupying  the  negative  X  axis,  and  the  optic  axis  in  the  medium  oriented  along  y. 


Fig.  10:  Coordinate  designation  for  the 
half-plane  diffraction  problem. 
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Since  <1  =  1,  8=0,  Eq.  (2-4)  reduces  to 


a2  a2  2 

— 1  +  €  — 2  ’  ko  €  1H^X*  y)  =  ■  iu,£n£  &  U-x')  6(y-y') 
bx  by 


The  boundary  conditions  are: 


|n  x  1£|  =  cos  6q  -  c  sin0o  =  0  on  S 


Radiation  condition  at  r  -  ® 
Edge  condition*  at  r  -  0 


(4-20a) 


(4-20b) 

(4-20c) 

(4-20d) 


The  physical  configuration  indicates  that  a  transformation  to  the  X-Y  coordinate 
system  will  simplify  the  task  of  solving  Eqs.  (4-20).  Applying  the  transformation 


^sin  v  -  cos  6  , 
o  o’ 


1  sin  6  /  \Yj 
o  o 


(4-21) 


to  Eqs,  (4-  20  a,  b)  yields: 


2  2  2 
+  J-^-T  +  2A.(S  )— 2— 

1  °  ax2  2  °  av2  3  0  axav 


+  k 


H( X,  Y)  =  L(2)  H  = 


i.ts  s  6(X-X')  6( Y- Y  ') 
o 


(4-22a) 


*The  edge  condition,  which  is  necessary  for  making  the  solution  unique,  is  a 
limit  on  the  allowable  singularities  of  the  field  components  near  the  edge.  The 
singularities  must  be  such  that  the  average  stored  electromagnetic  energy  in  a 
finite  volume  surrounding  the  edge  will  be  finite.  In  our  case,  the  energy  den- 

sity  is  W  =-r  E  .  -^(ojc).  E  +  p  |h|  .  The  components  of  E  are  derivable 

4L_  dx  -  -  o  J  2n  R  - 

from  H  via  Eqs.  (2-5),  thus  the  condition  W(r,  Qdrd  0<®  sets  an  upper 

&  S 

limit  on  the  growth  of  H  as  r  -  0.  In  the  process  of  solving  the  problem 
(Appendix  E)  this  is  a  crucial  point. 
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A3(0o)  ^c]h(X,Y)  =  l{1)H=0  on  j  ^»<X<  0 

(4-22b) 


The  radiation  and  edge  conditions  remain  unchanged,  because 


r  =l/777 


=  Vx2  +  Y2 


(4-23) 


In  Eqs.  (4-22) 


2  2 

Aj(w)  =  sin  w  +  £  cos  w 

A^(w)  =  cos^w  +  e  sin.^  w  =  £n(w)J^ 

A^(w)  =  (e  -  1)  sinw  cosw 


(4- 24a) 
(4- 24b) 
(4-24c) 
(4-24d) 


The  solution  of  Eqs.  (4-22)  proceeds  in  two  steps.  First,  plane  wave 
excitation  is  considered,  i.  e.  X'-®,  Y'-*m  and  X' /Y1  -  t,in0,(  From  the 
solution  of  that  problem,  the  solution  for  line  source  excitation  (finite  X'  and  Y  0 
may  be  synthesized.  In  both  cases  the  total  field  H  is  decomposed  into  incident 
and  scattered  parts  (IT  and  respectively) 


or 


H  =  H.  +  H 
1  s 


CP 


H 

i  U)  £  £ 
o 


cp.  +  :p 

1  s 


(4- 25  a) 
(4-25b) 


In  the  case  of  line  source  excitation,  cp.  is  a  cylindrical  wave  (Eqs. 
(3-13)  and  (4-11)).  By  letting  £,-<®  along  a  direction  B',  while  normalizing  the 
amplitude  to  unity,  the  cylindrical  wave  becomes  a  plane  wave  whose  ray  (energy 
transport)  direction  is  along  6': 

co  here  denotes  a  wave  (unction  and  should  not  be  confused  with  the 
same  symbol  employed  for  angles  in  other  Chapters. 
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with 


cp.  =  exp  £-iko(B2X  +  B3  Y) j 

B,  =  £e  sin  6q  sin  0/  +  cos  9q  cos  6  J 

B3  =  ^sin9ocos0,  -  e  cos  9^  sin  9'] 


(4-26) 

(4-26a) 

(4-26b) 


In  this  case,  cp  must  satisfy  Eq.  (4- 22a)  with  the  inhomogeneous  term  set  equal 
to  zero,  and  the  boundary  conditions  Eqs.  (4- 22b)  and  (4-20  c,  d).  From  Eqs. 
(4-22)  and  (4-25)  one  may  formulate  the  boundary  value  problem  for  the  scattered 


part  of  the  field  cp^: 


L{2)cp  =  0 
s 


L(i)  cp  =  -  cp.  on  S 

s  1 


(4-27a) 

(4-27b) 


By  using  Eqs.  (4-24)  and  (4-26),  Eq.  (4-27b)  may  be  written  explicitly  as: 

(1) 


L  \  =  iko[A3(8o)B2+A2(0o)B3]exPf-i  koB2  X]  = 

i k  s  sin(9  -S')  (  Y  n 

=  “JW^"  exP[-ikoB2X]  on  -»<X<0  (4‘27c) 


The  following  Fourier  transforms  are  now  defined: 

00 

*(a,  Y)  =  §  +  (a,  Y)  +  r_(a,  Y)  =  — |rcp(X,  Y)eiaX  dX 


2  n  _*> 


where 


Ma,  Y)  =  — Lr  f  cp(X,  Y)  eia  X  dX 

V2T  i 


(4-28a) 


(4-28b) 
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Mo-.  Y) 

©{>•:,  Y) 


Application  of  Eq.  (4-28a)  to  Eq.  (4-27a)  yields* 

d4* 


d“$ 


A. 


s 


2  dY 


Z*  2i“  A3lf  +  [*Z  -  Ml  *s  =  0  • 


After  assuming 


*  (a,  Y)  =  1(a)  e 

S 


i  q  Y 


(4-28c) 


(4-28d) 


(4-29) 


(4-30) 


substitution  into  Eq.  (4-29)  leads  to  the  quadratic  equation 


A  a2A  -k2 
ATq  +  - A— 

Ct  u 


0  , 


(4-31) 


the  solutions  of  which  yield  the  dispersion  relation 

7 


qi.2(a)  =  air  Wir(ko  -x 

2  2  V  c 


(4-32) 


e  ,  2 

The  function'y^-  |k  -  j  is  double  valued.  We  choose  that  value  for  which  the 

square  root  is  positive  when  real,  and  has  a  positive  imaginary  part  when  com¬ 
plex. 


*tn  the  following  discussion  A.(6  ),  A_(8  )  and  A_(6  )  will  be  written  simply 

I  O  CO  j  o 

as  Aj,  A^  and  A^  respectively.  The  argument  will  be  written  explicitly  only 

when  it  is  different  from  6  . 

o 
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The  corresponding  choice  of  q^  (upper  sign)  or  q2  (lower  sign)  is  discussed 
in  reference  (15)*.  In  order  that  cpg  satisfy  the  radiation  condition,  it  is 
shown  that  one  has  to  choose 


iq  Y 

ya,Y)=Ii(a)e  1 

for  Y>0  . 

(4-33a) 

iq2  Y 

*8(a,  Y)  =  I2(a)  e  c 

for  Y  <  0  * 

(4-33b) 

the  relation  between  Ij(a)  and  I2(a)  is  determined  from  the  fact  that  along 
Y  =  0 ,  the  tangential  component  of  the  electric  field  is  continuous:  i.  e.  bo*h 
L^cp  and  L^cpi  are  continuous  at  Y  =  0,  which  implies  that  L^'cp^  is  con¬ 
tinuous  along  Y  =  0.  Thus,  from  Eqs.  (4-33)  and  (4-28d)  one  gets 

[A2ql(°0  -  A3ajl1(a)  =  j^A£q2(a)  -  A3uJ  I2(a) 

which,  upon  use  of  Eq.  (4-3Z),  yields  the  relation 


Therefore 


Ij(a)  =  -  I2(a) 

iql  2{a)Y  < 

i  (a,  Y)  =  +  1(a)  e  c  for  Y$  0 
s  — 


(4-34) 


The  determination  of  1(a)  by  means  of  the  Wiener-Hopf  technique  pro¬ 
ceeds  as  in  reference  (4).  (Details  are  given  in  Appendix  E).  The  result  is 


1(a) 


lko 

Bzl 

2  - 

‘-Trrry 

(4-35) 


(a 


k  B, 

o  Z 


*  reader  is  cautioned  to  note  the  slightly  different  definitions  of  A.,  A_,  A, 
and  qj  2  in  reference  (15).  I  Z  3 
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from  which,  via  Eq.  (4-28d), 


<P3(X. 


--  aT 

r  b2  i 

*  da  exp  - 
j 

2  TT  ’O 

wyl 

(a- 

eT‘ko 


.  Y<  0  .  (4-36) 


Eq.  (4-36)  is  an  exact  solution  to  our  diffraction  problem.  It  should  be  noted 
that  it  has  been  derived  by  assuming  s  to  be  real  and  positive.  Next,  we 
investigate  the  behavior  of  the  integral  in  Eq.  (4-36)  when  e  is  complex. 

The  integrand  in  Eq.  (4-36)  has  branch  points  at  a  =  +  koN(9o)  which, 
for  real  and  positive  e,  lie  on  the  real  axis  of  the  complex  a  plane.  The 
branch  cuts  are  chosen  in  such  a  wav  that  on  the  entire  upper  Riemann  sheet  of 
the  double  valued  function^-^- /k^  -  ^-Y  the  imaginary  part  of  the  square  root  is 

_ _ _ . ^  i  • 


positive.  The  integration  path  shown  in  Fig.  11  lies  entirely  on  the  upper 
Piemann  sheet. 


Fig.  11:  The  path  of  integration  in  the  complex  a  plane 

a)  €  real  and  positive 

b)  complex,  !m  £>0 


47 


When  £  is  complex,  N{-  )  will  be  complex.  In  Chapter  III,  it  was  shown  that 
in  order  to  satisfy  the  radiation  condition,  N  has  to  be  positive  when  real  and 
ImN  positive  when  N  is  complex.  This  condition  is  satisfied  when 


0  <  arg  e  <  * 


(4-37) 


(see  Eq.  (3-30) ).  The  branch  points  then  move  into  the  first  and  third  quadrant 
in  Fig.  11(b),  and  will  not  interfere  with  the  integration  path.  Uniform  con¬ 
vergence  of  the  integral  in  Eq.  (4-36)  is  assured  if  as  ( 


Re  [i  q2  2(<Dy]  <0  ,  Y?Q  , 


(4-38) 


which  condition  is  seen  to  be  equivalent  to 


Im 


ijs  +_  A3^"’o^ 


A-A'-) 

c.  o 


>  0 


With 


£  =  |e| 


i  t 


(4-39) 


(4-40) 


and  using  Eqs.  (4-24),  condition  (4-39)  becomes 

VjeT  cos  vQ+  Is  1  sin^ro4_2^feTcos9o  sin0Q  sin-|y_>  0  \4~41) 

which  inequality  is  satisfied  for  any  v  and  £  satisfying  Eq.  (4-37).  One  may 
check  independently  that  the  individual  plane  waves,  and  therefore  the  total 
solution,  satisfy  the  power  radiation  condition  when  e<  0  (see  reference  15). 
Since  the  other  characteristics  of  the  solution  are  not  affected  by  letting  e  take 
on  the  range  of  values  menti'-  .ed  above,  Eq.  (4-36)  remains  valid  when  e  <  0. 

The  ray-optical  feature  of  this  result  may  now  be  determined  from  an 
asymptotic  evaluation,  the  details  of  which  are  given  in  Appendix  F,  for  positive 
or  negative  real  c.  The  method  of  steepest  descent  is  used,  and  it  is  shown 
that  there  arc  two  contributions  to  the  asymptotic  result,  one  from  a  saddle  point 
which  we  denote  by  v  and  one  from  a  pole  which  we  denote  by  <£>  : 

C  o 
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<P. 


=  cp  + 
c 


(4-42) 


It  is  also  shown  that  cp  exists  only  in  certain  regions  of  space  which  are 

S 

determined  by  the  relative  position  of  the  pole  and  the  saddle  point.  One  finds 
that  cp  is  that  part  of  the  solution  which  could  have  been  found  by  geometrical 
optics;  it  consists  of  the  geometrically  reflected  field  (obeying  the  reflection 
law  in  Eq.  (4-10))  and  another  pare  which  cancels  the  incident  field  in  the  region 
of  the  geometrical  shadow,  cp^  exists  throughout  space,  and  is  given  by 
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nk  rl4(6) 
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ik  rN(6)-i-£ 
o  4 


(4-43a) 


where 


f(ws) 


_  iy 
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esinG  sin9  +  cos9  cos9 
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esin9  sin9+cos6  cos6 

- N(9^(S)~? - 


c. sinG  sin0'  +  cos6  cos6’  e sin9  sinS+cosS  co«6 
N°GJ  NO7)  2  +  N(9j  N(6)  2 


(4-43b) 


For  e  =  1  (isotropic  medium)  this  formula  reduces  to  that  of 


Sommcrfeld^^ 


f(ws) 


s  =  1 


9-9'  5  -9 

2  sin - j sm — ^ — 

cos(C  -6')+  cost-  -9) 
o  o 


(4-43c) 


Eq.  (4-43a)  is  not  valid  in  the  neighborhood  of  the  poles  of  f(w  _).  Phys- 

s 

ically  these  poles  correspond  to  observation  point  locations  on  a  boundary 
between  geometrically  illuminated  and  shadow  regions  (see  Fig.  12).  J  or 
evaluating  the  field  in  these  transition  zones,  a  more  elaborate  calculation  of 
(4-36)  is  necessary^2<^  which  yields  a  result  in  terms  of  Fresnel  integrals.  It 
is  found  that  the  denominator  of  f(ws),  Eq.  (4-43b),  vanishes  at  two  angles: 


9 '  +  rr 


(4-44a) 
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9  2  =  tan'1| 


r  -( 1  -  s  tan^Q  )tanS  '  +  2  tan9 


2 - — 

(1-etan  9  )  +  2etan9  tan9 ' 
o  o 


(4- 44b)' 


A  comparison  of  Eqs.  (4-44b)  and  (4-10)  discloses  that  9^  is  precisely  the 
angle  of  the  limiting  reflected  ray  corresponding  to  rays  incident  from  the 
direction  3‘  (Fig.  12),  while  8j  describes  the  geometrical  shadow  boundary. 

It  has  therefore  been  confirmed  that  the  asymptotic  field  solution  in  the  exterior 
of  transition  regions  may  be  represented  in  simple  ray-optical  terms.  The  geo¬ 
metrical  optics  field  is  constructed  in  accord  with  the  discussion  in  Section  B 
of  this  chapter,  whereas  the  edge  gives  rise  to  a  diffracted  field  which  may  be 
interpreted  in  terms  of  diffracted  rays  (cf.  reference  3).  The  diffracted  rays 
progress  radially  outward  from  the  edge  and  therefore  have  a  radial  dependence 
as  given  by  the  factors  multiplying  f(w  )  in  Eq.  (4-43a)(see  also  Eq.  (4-16)  et. 
seq. ),  while  f(w  )  represents  the  "diffraction  coefficient"  which  yields  the 
starting  amplitude  of  the  ray.  t 


Fig.  12:  The  field  constituents  in  the  half-plane 
diffraction  problem. 
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The  plots  of  f(w  )  in  Figs.  13a,  b,  c,  d  show  the  influence  of  the  aniso- 
tropic  medium  on  the  diffraction  coefficient.  Whereas  the  diffraction  pattern  is 
symmetrical  in  the  isotropic  medium  (Fig.  13a),  the  presence  of  anisotropy 
may  introduce  substantial  distortion.  When  e  is  chosen  negative,  propa¬ 
gating  fields  are  confined  to  certain  angular  directions  so  that  the  diffracted 
field  emanating  from  the  edge  creates  its  own  region  of  illumination;  the 
pattern  function  in  the  shadow  zone  describes  the  angular  dependence  oi  the 
evanescent  diffracted  fields  ("evanescent  rays").  (It  is  to  be  kept  in  mind  that 
the  last  mentioned  regions  of  illumination  and  shadow  are  consequences  -of  the 
medium  properties  when  e  <  0  and  should  not  be  confused  with  illumination  or 
shadow  zones  caused  by  the  obstacle).  In  Fig.  13d,  the  incident  field  is  itself 
evanescent  but  gives  rise  to  a  propagating  diffracted  field.  This  aspect  is  of 
importance  in  connection  with  excitation  by  a  line  source,  to  be  considered  next. 

In  view  of  the  local  nature  of  high-frequency  (or  far  field)  diffraction,  it 
is  to  be  expected  that  the  amplitude  of  the  diffracted  field  is  proportional  to  that 
of  the  incident  field  at  the  edge.  Since  the  incident  radiation  is  locally  that  of  a 
plane  wave,  it  should  be  possible  to  construct  the  diffraction  field  due  to  arbi¬ 
trary  excitation  from  the  canonical  solution  (4-43a)  for  a  plane  wave.  This 

(3) 

argument  has  been  employed  in  the  geometrical  theory  of  diffraction  for 
electromagnetic  wave  propagation  in  isotropic  regions,  and  we  show  that  it 
remains  valid  also  in  the  anisotropic  case.  To  demonstrate  this  fact,  we  con¬ 
sider  the  case  of  a  cylindrical  wave  incident  on  the  half-plane,  the  rigorous 
solution  for  which  can  be  derived  from  that  of  plane  wave  diffraction.  The 
details  of  the  calculation  are  given  in  appendix  G.  We  assume  that  x7  and  y * 
in  Eq.  (4-20a)  are  finite  and  (x7/ y7)  =  tan07.  The  incident  field  cp.  is  now  in 
the  form  of  a  cylindrical  wave  as  in  Eq.  (4-11),  and  the  rigorous  solution  for 
the  secondary  field  cp^  is  shown  to  be  (see  Appendix  G): 

.i  ®  r-F-w',-  ifvv'1®  exp  -iaX+  iq,(a)Y 

+  1  „  exp  -i5  X  +  iq.  (c.)  Y  v  L  .  ^1  J 

cp  =  — dS - '•  1  - i  da — - -  . — & - 

8  8n  .i  ^k'oN(eo)-§  (*  +  $) 


(4-45) 
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for  Y'>0  and  Y<0.  After  carrying  out  an  asymptotic  evaluation  of  Eq.  (4-45) 

for  kQrN(0)»  1  (Appendix  F),  one  finds  again  that  the  solution  is  in  the  form 

of  Eq,  (4-42),  with  cp  representing  a  part  which  can  be  derived  by  geometrical 

8 

optics  (see  Appendix  C),  while  cp  -  the  diffracted  part-  is  given  by 
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(4-46) 


Since  f(w  )  is  the  same  coefficient  as  in  Eq.  (4-43b),  and  the  factor  distin- 
s 

guishing  Eq.  (4-46)  from  Eq.  (4-43a)  is  precisely  the  incident  ray  amplitude  at 
the  edge  (see  Eq.  (4-16)),  this  result  confirms  the  local  character  of  the  diffrac¬ 
tion  process  and  the  validity  of  the  geometrical  theory  of  diffraction.  The 
preceding  solution,  valid  for  positive  or  negative  e,  has  also  shown  that  the 
boundary  of  the  geometrical  shadow  (i.  e.  ,  the  boundary  of  the  domain  of  exist¬ 
ence  of  the  incident  wave)  is  given  by  the  ray  from  the  source  which  grazes  the 
edge,  regardless  of  whether  this  ray  is  propagating  or  evanescent.  Similarly, 
the  reflected  ray  boundary  is  the  one  predicted  from  geometric  optical  con¬ 
siderations  since  the  reflection  law  in  Eq.  (4-10)  remains  valid  even  when  both 
the  incident  and  reflected  rays  are  non-propagating  (see  reference  15);  in  the 
latter  case,  the  field  is  obtained  by  analytic  continuation  of  the  real  ray  solution 
to  imaginary  values  of  N(6).  The  Wiener-Hopf  representation  has  therefore 
furnished  a  simple  interpretation  of  the  "geometric-optical"  (better,  the  direct 
and  reflected)  field  when  e<  0,  and  it  has  shown  that  this  part  of  the  solution 
(and  its  spatial  domain  of  existence)  may  be  determined  from  direct  ray-optical 
considerations.  This  clarification  did  not  emerge  as  easily  from  the  alternative 
representation  mentioned  at  the  beginning  of  this  section.  The  diffraction  field 
in  Eq.  (4-46)  is  also  seen  to  be  in  a  form  which  is  readily  interpretable  for 
positive,  negative,  or  even  complex  z.  It  is  therefore  to  be  expected  that  the 
direct  and  reflected  field  contributions  in  the  presence  of  obstacles  of  more 
general  shape  may  be  constructed  by  the  techniques  of  geometrical  optics  even 

when  s  is  non-positive.  The  diffraction  field  for  a  variety  of  structures  may 

(2) 

be  determined  from  known  isotropic  results  by  Felsen1  s  method  when  z  is 
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positive,  and  analytic  continuation  should  then  provide  the  correct  behavior 
when  0<arge  <  tt.  These  conclusions,  valid  for  the  half  plane,  remain  to  be 
confirmed  for  other  configurations. 

Depending  on  the  location  of  the  half  plane  with  respect  to  the  .line  source 
a  number  of  interesting  situations  may  arise  when  e  <  0  so  that  ray  propaga¬ 
tion  is  limited  to  angles  6  such  that  |tan@|<  tan0c  =  wVuT.  Various  cases 
are  depicted  in  Figs.  14  (a-d).  In  Fig.  14  (a),  the  half  plane  is  in  the  shadow 
of  the  source  and  the  incident  field  near  the  surface  is  evanescent.  As  men¬ 
tioned  in  connection  with  Fig.  13,  the  edge  nevertheless  creates  its  zcne  of 
illumination  so  that  the  diffracted  field  in  this  case  dominates  the  evanescent 
reflected  field.  In  Fig.  14  (b),  a  portion  of  the  plane  is  illuminated  and  the 
edge  is  in  the  shadow.  The  reflected  field  is  now  the  same  as  from  an  infinite 
plane  and  the  diffraction  effect  is  small  since  it  is  caused  by  an  evanescent 
incident  wave.  When  the  edge  is  in  the  illuminated  region  of  the  source  but  the 
plane  is  mostly  confined  to  the  shadow  (Fig.  14  (c)),  both  the  reflected  and 
diffracted  fields  may  be  significant,  with  special  importance  assigned  to  the 
diffracted  wave  in  that  region  of  space  wherein  the  incident  and  reflected  waves 
are  evanescent.  Analogous  considerations  apply  to  the  final  case  wherein  the 
entire  plane  is  illuminated  (Fig.  14  (d)). 

D.  DIFFRACTION  BY  A  SMOOTH  CONVEX  CYLINDER:  POSTULATIVE 

APPROACH 

The  present  section  deals  with  diffraction  phenomena  caused  by  surface 
curvature.  It  is  known  from  studies  in  isotropic  media,  ..hat  such  effects  arise 
from  the  vicinity  of  the  shadow  boundary  on  the  scatterer,  and  the  diffracted 
fields  may  be  associated  with  waves  which  are  launched  at  the  shadow  boundary, 
travel  along  the  surface  into  the  shadow  region,  and  leak  energy  continuously 
during  their  progress  ("creeping  waves"!^2^"^^.  From  a  ray-optical  view¬ 
point,  the  creeping  waves  are  interpretable  in  terms  of  diffracted  rays  which 
are  excited  by  an  incident  ray  tangent  to  the  obstacle.  These  creeping  rays 
aiso  satisfy  a  generalized  form  of  Fermat's  principle,  and  the  mechanism  of 


Fig.  13:  The  diffraction  coefficient  f(wg) 


a)  e  =  l  (isotropic  medium) 

b)  e-  =0.  333 
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energy  transport  from  a  source  point  Q  to  a  point  P  in  the  shadow  region  may 
be  understood  from  the  ray  trajectories  sketched  in  Figure  i5.  For  a  ray  path 
QQj  P,  QQj  is  measured  along  the  incident  tangent  ray,  while  QjP^  is  the 
distance  traveled  by  the  diffracted  ray  along  the  surface  (with  an  associated 
leakage  of  energy).  Since  the  leakage  may  be  represented  by  rays  which  leave 
the  surface  tangentially,  the  path  segment  P^  P  is  the  associated  trajectory. 

If  the  sca'.terer  is  a  closed  cylindrical  structure  as  in  Fig.  15(b),  analogous 
considerations  apply  to  the  alternate  path  QQ^  P^  In  that  case  there  are 
additional  field  constituents  at  P,  arising  from  rays  which  have  encircled  the 
cylinder  one  or  more  times  before  shedding.  In  view  of  the  exponential  decay  of 
a  diffracted  ray  during  its  travel  along  the  surface,  these  latter  contributions  are 
frequently  negligible. 

The  preceding  features  have  been  incorporated  by  Keller^*^  into  a  pres¬ 
cription  which  allows  the  quantitative  construction  of  the  diffracted  field  in  iso¬ 
tropic  regions.  While  the  validity  of  these  postulates  has  not  yet  been  established 
in  general,  the  successful  verification  for  various  special  cases^^^S)  jeR<jg 
strong  support  for  their  applicability  to  smooth  convex  surfaces  of  relatively 
arbitrary  shape.  We  shall  therefore  follow  similar  arguments  to  obtain  repre¬ 
sentations  for  the  diffracted  field  when  the  scatterer  is  embedded  in  an  aniso¬ 
tropic  medium.  Referring  again  to  Fig.  15,  the  incident  field  at  Qj  caused  by 
a  line  source  at  Q  is  given  by: 


where  A^  is  the  reference  amplitude 


-L  i 

"k  n(s:) 

o  i 


(4-47a) 


We  denote  the  field  amplitude  associated  with  the  diffracted  ray  near  the  cylinder 

surface  by  Ad(0.  where  >  is  the  arc  length  measured  from  Qj.  To  relate 

A,  to  A.,  assume 
d  l 
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Ad(0)  =  IXQ^A.tQj)  , 


(4-48) 


where  D(Oj)  is  a  diffraction  coefficient  which  is  a  function  of  the  properties  of 
the  surface  and  the  medium  at  Qj.  Because  of  the  continuous  leakage  of  energy 
from  the  surface,  it  will  be  further  assumed  that 


AdU)  =  Ad(0)  exp 


* 

-f  ad'  , 

*  o 


(4-49) 


where  a  is  a  decay  coefficient  which  also  is  a  function  of  the  local  properties 
of  the  surface  and  the  medium.  If  A(d)  denotes  the  field  amplitude  along  the 
ray  which  leaves  the  surface  tangenti?iiy  at  Pj,  then  we  assume  that 

A(0)  =  D(P j)  Ad(0)  exp  -J  1  a  d  A  ,  (4-50) 

Lq.  J 


where  D(P^)  is  again  a  diffraction  coefficient.  As  long  as  the  field  satisfies 
conventional  reciprocity  conditions  (which  is  true  for  the  uniaxial  case),  we  may 
assume  that  D(Qj)  and  D(Pj)  have  the  same  dependence  on  the  surface  and 
medium  properties  at  and  Pj ,  respectively.  In  cases  of  more  general 
anisotropy,  this  point  has  to  be  reexamined,  and  more  general  assumptions  have 
to  be  made  in  accord  with  the  reciprocity  relations  satisfied  by  the  medium. 
Also, 
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(4-51a) 


The  field  at  P  due  to  one  ray  is  therefore  given  by 

.  . /  TP. 
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where  Y(P)  is  the  phase  change  along  the  trajectory: 

Pj 

Y(P)  =  ikQ  jdj  N(6J)  +  djNtOj)  +  J  N[9U)]dt 


(4-52a) 


By  analogy  to  the  isotropic  case,  we  may  expect  the  possible  existence  of 
different  "modes"  having  different  decay  exponents  and  diffraction  coeffi¬ 
cient  0^,  p  =  i,  2,  .. .  .  The  total  field  at  the  point  P  in  Fig.  15(a)  is  there¬ 

fore  assumed  to  be  given  by  a  sum  of  all  of  these  modes, 
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(4-53) 


If  the  scatterer  has  a  closed  cross  section  (Fig.  15(b) ),  one  has  to  add  the  con¬ 
tributions  from  rays  which  have  encircled  the  cylinder  m-times,  m  =  l,  2 . 

The  formulas  for  the  corresponding  rays  differ  from  that  in  Eqs.  (4-52)  in  that 
the  integration  interval  extends  from  Qj  to  (P^  +  m  L)  and  from  to 
(P^+  mL),  where  L  is  the  circumference  of  the  cylinder.  Using  the  formula 
for  summation  of  an  infinite  geometric  progression,  we  can  write  the  sum  of  all 
these  rays  in  closed  form: 


(4-54) 


Thus,  the  formula  analogous  to  Eq.  (4-53),  which  gives  the  total  diffracted 
field  in  the  shadow  region  of  a  closed  cylinder,  is: 
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In  these  detailed  formulas  (Eqs.  (4-53)  and  (4-55) ),  it  is  often  sufficient 
to  retain  only  the  p=l  mode  (having  the  lowest  decay  exponent  d^)  and  to 
ignore  the  denominator  term  in  the  summand  (i.  e.  to  neglect  the  rays  which 
have  encircled  the  cylinder  completely  one  or  more  times).  This  results  in  an 
expression  comprising  one  diffracted  ray  for  the  field  in  the  shadow  of  an  open 
cylinder  (as  in  Eq.  (4-52) ),  and  two  diffracted  rays  for  the  field  in  the  shadow 
of  a  closed  cylinder. 

While  the  preceding  discussion  constitutes  in  effect  an  application  of  the 
geometrical  theory  of  diffraction^^ to  an  anisotropic  medium,  and  is  there¬ 
fore  based  on  a  series  of  postulates  similar  to  those  in  the  isotropic  case,  con¬ 
firmation  of  its  validity  may  be  obtained  directly  for  special  configurations,  when 
the  medium' is  uniaxial.  When  8  =0,  with  a  and  C  real  and  positive,  the  scat¬ 
tering  problem  in  the  uniaxially  anisotropic  region  is  reducible  to  an  equivalent 
isotropic  problem  by  the  previously  described  method  of  coordinate  scaling 
(Eq.  (2-6)).  Since  the  isotropic  analogue  of  Eq.  (4-53)  and  (4-55)  has  been  con¬ 
firmed  for  such  obstacles  as  circular,  elliptic  and  parabolic  cylinders^ 
an  application  of  the  scale  transformation  yields  a  corresponding  asymptotic 
solution  for  the  uniaxial  case.  It  may  be  verified  that  the  result  derived  in  this 
manner  furnishes  not  only  the  general  form  given  in  Eqs.  (4-53)  and  (4-55)  but 
also  yields  specific  expressions  for  the  coefficients  D  and  &  . 


To  illustrate  these  remarks,  consider  the  equation  for  a  single  ray  vEq. 
(4-52))  in  an  isotropic  medium  (uv  -  coordinate  space)  with  wave  number 

k  =  k  Ja  =  k  lax  : 

o’  o^ 


r  - 

P, 

_  _  i»  i _ 

P  1  _ 

ik 

d  +  d  +  I  d<. 

-  J  apdt 

^1  J 

d' 

and  df.  defined 

in  Fig.  16. 

(4-56 


Fig.  16:  The  uv  -  coordinate  space. 


wtutmmr 
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D  and  a  in  the  isotropic  medium  have  been  calculated  by  Keller 

.n/12  1/6  1/6 

_  n  e1  kQ  [b(u.v)] 

p  =  2  a1/12  ,1/4  7T75 


(16). 


'[W] 


(4-56a) 


a  =-ik  1/3^/6 
P  o 


[b{u,  v)] 


•2/3 


(4-56b) 


where  b(u,  v)  is  the  local  radius  of  curvature  of  the  perfectly  conducting 

scatterer,  and  the  numbers  7  are  defined  by 

P 

in/3  ,-1/3  .  ,  . 

T  =  q  e  6  (4-56c) 

P  P 


with  the  numbers 


denoting  the  roots  of  the  equation 


^A{q)"  J  cos(w3  -qw)dw=0 

o 


(4-56d) 


By  applying  the  transformation  x=q^X,  y =v/e\  the  resulting  solution  in  the 
x-y  space  is  known  to  satisfy  Eq.  (2-4),  the  boundary  condition  E^O  on 
the  trasformed  obstacle,  and  the  radiation  conditon  at  infinity.  Since 


ii I  =JE  =Jz 

du  dx  ’€ 


cot  5, 


where  5  is  measured  from  the  positive  y-axis,  one  obtains  for  the  transformed 
length  element 


d-f  (u,  v)  -<£77^  =71  dx^ +  dy 


1 1  .  2  1  2 

T<lzdx  +  ed  y 

Vc  2  ■ 2 


fdx  +  dy 


=  —  N(6)  dt(x,  y) 

$ 


(4-57a) 
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where  N(0)  is  the  ray  refractive  index  (see  Eq.  (4-8))  and  A=det|e|=  ae 
(since  8=0).  In  the  same  way  one  can  show  that 


kd(u,  v)  =  koN(9)d(x,  y) 


kd'(u,v)=  kQ  N(9  )d  (x,  y)  . 


(4-57b) 

(4-57c) 


The  phase  function  in  Eq.  (4-56)  is  seen  to  transform  into  the  one  given  in  Eq. 
(4-53),  with  due  cognizance  taken  of  the  fac.t  the  N(9)  varies  over  the  part  of 


the  trajectory  on  the  obstacle  surface.  Similarly. 


[b(u,  v)] 


■  2/ 3  \du 


N^(S) 


(4-58) 


where  b(x,  y)  is  the  radius  of  curvature  of  the  obstacle  in  the  x-y  space.  The 
decay  exponent  in  Eq.  (4-56)  is  then  converted  into: 


1 


f  a  dt  =  -  i 
J  p 


1/3  . 1/3 


ik  "'A* 
o 


r  f  dj 

Po  RiiF^N^U) 


from  which  one  may  identify 

i  k  1/3  L1/3  - 


a  =  - 
P 


t — m — pi 

b{l)\cii  N[8R) 


(4-59) 


The  diffraction  coefficient  is  obtained  from  in  Eq.  (4-56a)  upcn 
inserting  from  Eq.  (4-58). 

Thus  it  has  been  shown  that  the  field  postulated  in  Eqs.  (4-53)  or  (4-55), 
together  with  the  results  for  Ct^  and  from  above,  describes  correctly  the 


f 
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diffracted  rays  on  simple  convex  surfaces  embedded  in  a  uniaxially  anisotropic 
medium  characterized  by  real  and  positive  a  and  e.  In  the  derivation  of  the 
result,  it  was  assumed  that  the  asymptotic  field  representation  in  the  aniso¬ 
tropic  medium  may  be  derived  by  applying  the  scale  transformation  to  the 
asymptotic  field  solution  in  the  equivalent  isotropic  configuration.  The  validity 
of  this  procedure,  not  really  in  question  when  e>0,  a>0,  is  confirmed  for  two 
special  cases  treated  in  the  next  sections. 

The  series  in  Eqs.  (4-53)  and  (4-55)  remain  convergent  when 
0<arge<Tt ,  as  may  be  verified  by  examining  the  exponents  in  Eqs.  (4-53)  and 
(4-55),  which  retain  a  negative  real  part  that  increases  with  p  for 
0<  argN(8)<  n/2  (corresponding  to  0<arg£<n).  It  is  plausible,  therefore,  to 
expect  the  preceding  formulas  to  hold  in  this  extended  parameter  range.  A  con¬ 
firmation  of  this  statement  does  not  follow  immediately  from  the  previously 
utilized  scaling  technique  since  the  scaled  coordinates  are  non- real  when  e  is 
non-positive.  By  examining  the  problem  of  diffraction  by  a  parabolic  cylinder, 
however,  the  solution  may  be  phrased  in  a  form  which  validates  the  analytic 
continuation  of  the  asymptotic  formulas  to  non-real  e.  The  analysis,  which 
establishes  the  validity  of  Eqs.  (4-53)  and  (4-55)  for  0<arg€<n,  is  presented 
in  the  next  Section. 


E.  DIFFRACTION  BY  A  PARABOLIC  CYLINDER:  RIGOROUS  ANALYSIS 


In  order  to  confirm  the  general  results  which  were  derived  in  the 
previous  section  in  a  postulative  way,  the  rigorous  solution  to  a  special  case 
is  worked  out  in  this  section.  Felsen's  method^  is  used  to  obtain  the  expres¬ 
sion  for  the  field  in  presence  of  a  perfectly  conducting  parabolic  cylinder 
embedded  in  a  uniaxially  anisotropic  medium  (8=0,  (X=i  in  Eqs.  (<I-1)).  The 
parabolic  cylinder  surface  S(x,  y)  is  given  by  the  equation 


Y  = 


(4-60) 


with  the  optic  axis  parallel  to  the  y  axis.  This  obstacle  has  a  varying  radius 
of  curvature,  with  the  smallest  radius  of  curvature  b  at  x=0.  Also,  the  angle 
$  between  the  surface  normal  and  the  optic  axis  varies  continuously.  A  mag- 


■11— 
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netic  line  current  is  located  at  the  point  Q  to  the  left  of  the  obstacle  (x<0). 

As  shown  in  Chapter  II,  Section  A,  the  magnetic  field  has  only  a  z-component 

H(x,  y)  from  which  the  two  electric  field  components  E  and  E  may  be 

x  y 

derived.  H(x,  y)  is  defined  uniquely  by  the  equations 


'  *2  .2  \ 

_!_+  e-^  +  k^  H(x,  y) 

9y  j 


-i'X'Co€  6(x-x/)  b(y-y')  , 


(4-6la) 


on  S 


(4-6lb) 
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Radiation  condition  at  rH®  ,  (4-6lc) 

with  k=ko/e" ,  and  $=tan  *  ^  denoting  the  angle  between  the  normal  n  and 
the  optic  axis. 

It  is  convenient  to  define 

y  =  y+a  ,  (4-62a) 

where  a  is  an  arbitrary  (possibly  complex)  finite  constant.  The  equation  of 
*  * 

the  surface  S(x,  y)  becomes 

2 

y  =  a-fr-  .  (4-62b) 


It  is  seen  that  the  change  of  variable  (4-62a)  does  not  affect  the  quantities 


■1  dy 
dx 


so  that  in  Eqs.  (4-6la)  and  (4-6lb),  y  can 


^  ,My-y')  and  <|f=tan 

y  3y 

A 

be  replaced  by  y  and  S  by  S  .  The  radiation  condition  requires  outward  flow 

n — ? 

of  energy  as  r=1fx  +y  -  •  .  As  "a"  is  a  finite  quantity,  the  radiation  condi¬ 
tion  Eq.  (4-6lc)  also  remains  unchanged. 

If  we  choqse 

a  =  €/lb  ,  (4-63) 


then  upon  introduction  of  the  transformation 

x  =  u,  y  =  <fe  v 


the  equation  for  S  becomes 


=  _ 

2  2  b  ' 


(4-64a) 


(4-64b) 


which  is  the  equation  of  a  constant  coordinate  surface  in  a  parabolic  coordinate 
system.  Via  Eq.  (4-64a)  the  boundary  value  problem  in  the  u-v  coordinate 
system  reads 


G(R.R')  =  -6  (u-u')  6(v-v') 


(4-65a) 
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dG  _  3G  .  l  3G  l  .  Je  b 

— x-  =  3—  sin  y  --5 —  cos  v  =  0  on  v  =  --—■■■ 
Sn  ou  3v  2 

Radiation  condition  at  R 


u 


2  Jt  b 


=4?+v*  -•  , 


where 

and 


G(R.R')  = 


H 


♦  £  •  *“-*£& 


A  set  of  parabolic  cylinder  coordinates  is  now  introduced  via 

U  =  T)  §  . 

1  In1  ~2\ 

v  =  j  (ri  -  s  )  , 

the  range  of  variation  of  t]  and  S  being 

0  <T)  <  +  00 
-«°<S  <+  ® 


(4-65fc) 

(4-65c) 

(4-65d) 

(4-65e) 

(4-66a) 

(4-66b) 


(4-66c) 

(4-66d) 


iVe  assume  e  to  be  real  and  positive  at  this  stage.  Thus  y,  u,  v,q  and  §  are 
ill  real.  Later  on  we  will  let  e  be  complex  and  investigate  analytic  continua- 
;ion  of  G  in  the  complex  e  plane.  In  the  T}-s  space, Eqs.  (4-65)  become 


v  +  4+kv+s2) 

dn  SI 


G  =  -6(r,-T!')  6(§-§  ')  , 


=  0  on  rj  =  r|o  =  b  , 


12-2 

Radiation  condition  at  R=-^  fa  +  5  )  ~® 


(4-67a) 

(4-67b) 

(4-67c) 


The  surface  ri  =  Tl0  is  real  only  when  €  is  real  and  positive.  For  complex  or 
negative  real  e  it  cannot  be  drawn  as  in  Fig.  18,  but  Eqs.  (4-67)  remain  valid. 
The  surface  has  a  large  radius  of  curvature  if  b  is  large  with  respect  to  wave¬ 
length.  T)  is  related  to  b  via  Eq.  (4-67b),  so  we  will  assume 
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rj  >  -  £  >  0  1}  >  £  >0 


Fig.  18  Parabolic  cylinder  coordinates. 

|kr^2j»l  ,  (4-68) 

with  the  absolute  value  sign  in  Eq.  (4-68)  applying  when  €  is  not  real  and 
positive.  A  representation  of  the  Green's  function  G  in  terms  of  parabolic 
cylinder  functions  will  now  be  found  by  means  of  the  characteristic  Green's 
function  method1  .  According  to  that  method  G  is  given  by 

G(R,R')  =  §  g^n.^';  X)g§(5.5,;\)  dX  (4-69) 

where  X  is  the  separation  constant  of  the  partial  differential  operator  in  Eq. 

(4-67a),  and  the  functions  g  and  gr  are  solutions  of  the  one -dimensional 

*»  ^ 
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equations 


d  +  .2J2  1 

— j  +k  T)  4 

dn 


-^y+kV+X 

d5  L 


X)  =  -6(ti-ti')  , 


=  -6  (§-§') 


(4-70a) 


(4-70b) 


The  contour  of  integration  in  Eq.  (4-69)  closes  at  infinity  and  surrounds  all  of 
the  singularities  of  either  g^  or  g§  in  the  complex  X  plane.  The  boundary 
conditions  ok  g^  and  g,  are  identical  with  those  satisfied  by  G  in  the  q  and 

9  domains,  respectively.  Two  linearly  independent  solutions  of  the  homogeneous 
equation  (4-70a)  are  given  by 


fjOilX)  =  D  l  n  (qVTTk)  . 
"2  ‘2k 


(4-71a) 


f2(ri;X)  =  D  x  .^C-q/TTk)  , 
"2  "2k 


w 


ith  fT  =  e"iT™ 


(4-71b) 


where  the  parabolic  cylinder  function  Dv(z)  satisfies  the  differential  equation 

(4-71c) 


.2  i  2 


Dv(z)  =  O 


To  avoid  ambiguity  in  the  asymptotic  expressions  of  fj  and  f  for  |x|-»  , 
it  is  necessary  to  define  the  relations  between  arg  v  and  arg  X 


»rg  v  =  argX-1.  when  |v|-«  and  argv<0 


arg  v  -  arg  X  +  ■—  when  |v|-»  and  arg  v>  0 


(4-71d) 
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Also,  the  Wronskian  is  given  by 


wfi-‘4  =  vi-vi 


2  j-ink 


(4-71e) 


Only  the  function  fj  satisfies  the  radiation  condition  at  rj-**.  This  maybe 
seen  from  the  asymptotic  formulas  for  large  values  of  |z|  which  are 


Dv(2)~zVe'7'2/4  [l+OU'2)]  , 


|  z|»|v|  and!  arg(z)j<2il  , 


(4-72a) 


Dv(z)~zVe-z2/4[l  +  0(z-2)]  z-v*lcinv^/4[i+0(l'2,] 


Ui^lvj  and  -^>arg(z)>^  . 

Since  t]  is  real  and  positive,  one  finds  from  Eq.  (4-72a) 

.1  _AA  z 

D  l  iX(rltfrnk)~(7ltf2ik)  2  2kelkTl/ 2  , 


(4-72b) 


(4-72c) 


2  2k 


ik  R 


which  is  proportional  to  e  as  r]-*®  .  In  view  of  Eq.  (4-72b),  the  same 
relation  is-not  obtained  for  f^.  A  linear  combination  of  f^  and  f 2  which 
satisfies  condition  (4-67b)  is  given  by 

D'  .  ..(-rih) 

-2  1 x  O  T 
2  2k 


1 


f3(-:;U  •  D 

"  2  "  2k  j 


0  *  n(”<h,*T-rn^ 


2  2k 


I  "Ik 


with 


and 


h  =  V-2 ik 
3 


Dv(axo)S  ^DJ3x)|x-x 


(4-73) 

(4-73a) 

(4-73b) 
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•jn 

--2~<arg  X  < -tt  we  have  0<  arg  v<  —  .  and  the  lower  sign  in  Eq.  (4-7 ha.)  is 
chosen.  In  this  way  Eqs.  (4-74a,  b)  yield: 


"2Vj  , -n<arg  X<2  (4-76a) 

^ln-n'l.e  VUn-V-znj]  .  w<arg  x<.„  (4.76b) 

which,  in  view  of  the  fact  that  Tj+ri^Zn  ,  is  exponentially  small  throughout  the 

0  TT 

indicated  region.  Special  care  is  needed  in  the  interval  -TT<arg  X<  -  -=■  which 

3n  ^ 

corresponds  to  --^-<arg  v<  -tt.  In  that  range  Eq.  (4-75b)  cannot  be  used 
directly.  Nevertheless  one  can  define  v  '  =  |  v'  |  e*^  ,  where  =  arg  v  and 

find  the  asymptotic  behavior  of  T(v)  from  the  definition  of  that  function^*^. 

By  doing  so,  Eq.  (4-76a)  is  shown  to  hold  over  the  whole  region  of 

-n<arg  X<  ^  .  Similarly,  two  independent  solutions  of  the  homogeneous  equation 

(4-70b)  are  given  by 


hA(S ;\)  =  D  j  u(W-2ik) 
~2+2k 


(4-77a) 


h2(5;X)  =  D  1  u(-§  V-2ik) 

"  2+2k 


with 


W[hi,h2]  =  h1h2-  h2hL  -  jf{  ITS 

\Z  '  2k/ 


(4-77b) 


(4-77c) 


i  k  R 


From  Eq.  (4- 72a)  it  is  easy  to  see  that  h.  is  proportional  to  e  as 
§-  +  °°,  and  h2  is  proportional  to  as  •a-*-*’  ,  thus^^ 


h ,(L;  \)hJ%  :X) 

wjh^— 


(4-77d) 

with  §>  and  defined  as  in  Eqs.  (4-74).  The  exponential  decay  of  g^  as 
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I  X  |-*«  can  be  shown  in  a  way  similar  to  the  one  used  for  g^  . 

The  singularities  of  in  the  complex  X  plane  arise  from  the  poles 
of  the  f  function,  has  a  series  of  poles  at 

=  -ik(2n+l)  ,  n  =  0,l,2... 

does  not  have  poles  at  =  i  k(2n  4  1),  because  at  these  points 


D'i 
2  "  2k 


D  1  U('T1<h)+D  1  i\(T,<h)  15T 

"I  '  2k  "  ’  ’  1  xK 


2  2k 


■RTF 


2  2k 


=  0 

x=\ 

n 


This  can  be  shown  by  using  the  relations 

D>>  =|Dv(Z)-Dv  +  1(z)  , 

valid  for  any  complex  v  and  z,  and 


Dv(-z)  =  (l)VDv{z),  v  =  0.1.2.. 


Nevertheless,  g^  has  poles  in  the  complex  X  plane  at  those  points  where 


D  .  iX  (n  h)  =  0 
i  p  o 

'2  "2k 


(4-78a) 


The  values  X^  have  been  calculated  by  Keller  and  Levy 

totic  values  for  kr,2»  l  are 

o 

X  ~kV  [l+2T  (kri2)'2/3] 

P  o[  pv  o'  J 


(32) 


Their  asymp* 


(4- 78b) 


where  is  defined  by  Eqs.  (4-56c,  d).  The  path  of  integration  runs  along  the 
real  axis  from  X=-®+is  to  X=  +  *+is  with  -k<  s  <  0 ,  and  closes  by  a  semi¬ 
circle  on  either  the  left  or  the  right.  Due  to  the  exponential  decay  of  and 
g?  as  |x|-»®  ,  the  semicircle  does  not  contribute  to  the  integral.  Defining 
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U  =  -  — +  — 

2  2k 


(4-79a) 


and  using  the  relation 


n-v-mi  +  v)  ~ 


sinn  v 


(4- 79b) 


Eq.  (4-69)  yields  explicitly 


G  =  G1  +  G2  = 


-  — +  i® 
2 


*i  I  7ii^r  %tlV-i<t|D-ii-i|sh|k-i,-1<t|t D-u - 1 (T1<h) g;~7(Tr£r 

1  ~r*  **  *  O 

_  A  .  X«  *- 

2  (4-80) 
Eq.  (4-80)  is  a  rigorous  solution  of  the  problem  formulated  in  Eqs.  (4-65)*. 

The  solution  of  the  original  anisotropic  problem  stated  in  Eqs.  (4-6l)  is  obtained 
by  expressing  the  variables  in  terms  of  x  and  y.  From  Eqs.  (4-66), 

(4-64),  (4-63)  and  (4-62)  one  finds 


’"•JIT  =  */k(R+v)  =  "\/kl/x“+  y  2/e  +  y/e  = 


ik«  h»y2+ex^+'i?+“^7  f  y  +  4i 


(4-81a) 


ijk  =  Vk(R-v)  Wk 


5  b 


€  C 

2  'y_  2b 


(4-81b) 


If  the  square  roots  are  defined  as  positive  when  real  and  having  positive  imagi¬ 
nary  parts,  then  it  can  be  shown  that  for  0<arg£<~ 


-^lar8  fa  -4^2  i  k) 


(4-82a) 


n  .  —  *  -  ,  I  .  n 

-  4larg('-  V-2  ik)<- 


(4- 82b) 


Rice'  and  Jones'  ;  have  discussed  the  diffraction  of  plane  waves  by  a 
parabolic  cylinder.  Eq.  (4-80)  can  be  shown  to  agree  with  their  expressions 
when  iR,|-®.  They  did  not  use  the  characteristic  Green's  function  method 
for  deriving  their  results. 
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It  can  be  verified  by  inspection  that  Eq.  (4-82a)  follows  from  Eq.  (4-8la).  The 
derivation  of  Eg.  (4-82b)  from  Eq.  (4-81b)  is  less  obvious.  One  has  to  show 
/  2  2  v  £  ^  t 

that  Imyy  +  €x  +  e  b +  4b  —  2b  *  This  can  be  done  by  observing  the  fact  the 
x  and  y  are  confined  to  the  region  exterior  to  the  parabolic  cylinder. 

Up  to  this  point  e  has  been  taken  as  positive  real,  thereby  implying 
positive  real  iq,  rj '  and  rjQ.  For  a  source  location  as  in  Fig.  18,  %'  is  negative 
real,  whereas  §  may  be  positive  or  negative  real.  Also  arg  (ti</-2  ik  )=  -  , 

arg  (-T1 V-2  ik)  =~  ,  and  similarly  for  the  5  variable.  For  bounded 
v  and  |x|~*®,  D^(x)  is  exponentially  small  in  |arg(x)|<-^  and  exponentially 
large  in  j  arg  (-x)  |  <^  .  It  is  seen  from  Eqs.  (4-82)  that  if  e  is  complex  and 
Im  £>0,  the  asymptotic  behavior  of  D^(x)  does  not  change,  and  the  considera¬ 
tions  concerning  the  radiation  condition,  which  lead  to  the  solution  in 
Eq.  (4-80),  remain  valid.  The  convergence  of  the  integral  (4-80)  can  also  be 
shown,  by  using  Eqs.  (4-82)  and  (4-75).  Thus,  since  the  solution  (4-80)  of  Eqs. 
(4-67)  remains  equivalent  to  that  of  Eqs.  (4-6l)  when  e  is  complex,  use  of 
Eqs.  (4-62),  (4-63)  and  (4-64)  permits  the  construction  of  the  solution  in  the 
anisotropic  medium  for  Im  C>0. 

The  asymptotic  expansion  of  Eq.  (4-80)  is  different  in  the  various  phys¬ 
ical  regions  of  figure  19.  If  the  observation  point  lies  in  the  illuminated  region, 
a  direct  ray  from  the  source  as  well  as  a  ray  reflected  from  the  obstacle 
according  to  the  laws  of  geometrical  optics  contribute  to  the  field.  The  direct 
ray  is  given  by  in  Eq.  (4-80).  The  reflected  ray  may  be  derived  by 
applying  the  saddle  point  method  of  evaluation  to  the  integral  G^  in  Eq.  (4-80). 
If  the  observation  point  lies  in  the  shadow  region,  the  first  order  asymptotic 
evaluation  of  G^  yields  an  expression  similar  to  the  asymptotic  expression  of 
Gj,  but  with  an  opposite  algebraic  sign.  This  indicates  the  cancellation  of  the 
field  to  the  first  order.  A  more  detailed  calculation  of  G  in  terms  of  the 
residues  in  the  complex  U  (or  X)  plane  lends  itself  to  an  interpretation  in  terms 
of  diffracted  rays.  If  the  observation  point  lies  in  the  transition  region  near  the 
boundary  between  the  illuminated  and  the  shadow  zones,  both  methods  fail. 
Actually,  in  this  region,  the  field  cannot  be  expressed  in  simple  ray  terms.  A 
numerical  evaluation  of  the  integrals  yields  a  smooth  transition  function  which 
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Fig.  19  Diffraction  by  a  parabolic  cylinder: 

The  illuminated,  transition  and  shadow  zones. 

connects  the  different  asymptotic  forms  on  both  sides  of  the  transition  so:iC 
because  the  shadow  boundary  does  not  introduce  any  physical  discontinuity  into 
the  fields. 

The  asymptotic  expression  for  the  diffracted  field  in  the  shadow  zone  is 
calculated  as  follows: 

The  only  poles  of  G  to  the  left  of  the  path  of  integration  in  the  complex  Ji  plane 
are  given  by  Fq.  (4-78).  The  path  can  be  closed  by  an  infinite  semicircle  en¬ 
closing  the  left  half  plane,  and  the  integral  is  given  exactly  by  the  sum  of  resi¬ 


dues: 
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(4-83) 


If  the  source  point  is  removed  to  infinity  and  the  result  is  normalized  in  such 
a  way  that 

H^(k|_R-R/ 1 )  -*exp  [-  ik(u  sinip  '  +v  coscp  ')]  ,  then  by  using  Eqs.  (4-66)  and 
(4-72)  it  can  be  shown  that 


v-h§  >rVi(hT1  > 


i  2^2n  cos^5  / 2 


/  u 


;  tan  cp  =  — 


The  resulting  Eq.  (4-83)  agrees  with  the  expression  derived  by  Rice^^  for  the 
case  of  plane  wave  diffraction. 

Keller^^  has  shown  that  the  exact  expression  given  by  Eq.  (4-83)  can 
be  written  asymptotically  as 


i 

exP  [ik^'+ d+  J  d*.)] 


•  l  i5P<F,1yQ)  exp  [iTPkl/3  jcw2/3<4 

p=o  P  J 


(4-84) 


with  b  denoting  the  radius  of  curvature  along  the  diffracted  ray  path,  and  ^  • 

p 

the  diffraction  coefficient  given  by  Eq.  (4-56a)  (see  Eq.  (4-56)  and  Fig.  16). 

It  is  noted  that  for  complex  values  of  e,  these  quantities  (as  well  as  r\,r\' 
and  §’  which  are  related  to  them)  no  longer  remain  real.  Nevertheless,  the 
calculation  leading  from  Eq.  (4-83)  to  Eq.  (4-84)  remains  valid,  because  (as 
shown  in  the  discussion  following  Eq.  (4-8l|  the  same  asymptotic  forms  of  the 
parabolic  cylinder  functions  may  still  be  used  when  their  arguments  become 
complex  as  long  as  Im  e>0.  One  may  now  write  the  asymptotic  expression 
for  G,  Eq.  (4-84),  in  terms  of  the  original  x-y  coordinates  by  changing  back 
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via  Eqs.  (4-81).  The  expression  so  obtained  is  exactly  the  one  predicted  in  the 
preceding  section  (Eq.  (4-53^,  and  therefore  justifies  use  of  the  formula  for 
complex  6 .  It  is  also  shown  in  reference  (16)  that  the  series  in  Eq.  (4-84)  is 
equivalent  to  that  in  Eq.  (4-83)  only  if  the  observation  point  lies  in  the  shadow 
region.  In  the  illuminated  region,  one  performs  a  first  order  saddle  point 
evaluation  of  the  field  using  directly  the  integral  representation  of  Eq.  (4-80). 
For  finding  the  saddle  point  locations,  it  is  most  convenient  to  substitute  into 
Eq.  (4-80)  the  asymptotic  forms  of  the  parabolic  cylinder  functions  which  are 
obtainable  by  applying  the  WKBJ  method  to  Eqs.  (4-70).  The  results  are: 


■  x(ti  V-2ik)~ 
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1/4  exp 
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(4-85a) 
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(4-85b) 


2  /  2  2  /  2 

These  solutions  reduce  to  Eq.  (4-72)  if  rj  »X/k  and  §  =5»X/k  ,  respectively, 

and  must  therefore  have  at  least  the  same  region  of  validity  in  the  complex  planes 

o{{t\4-2  ik)  and  (§  V-2ik).  This  assures  us  that  the  asymptotic  evaluation  of 

the  field  can  be  performed  in  the  u-v  space,  and  the  result  then  expressed  in 

(24) 

terms  of  x;  y.  The  calculation  has  been  carried  out  by  Jones  for  the  case 
of  plane  wave  diffraction.  Our  case  differs  only  in  a  few  details,  which,  are 
discussed  in  Appendix  H. 


F.  RIGOROUS  ANALYSIS  FOR  ANOTHER  SPECIAL  CASE:  A  CLOSED 
CYLINDER 


Diffraction  by  a  large  elliptic  cylinder  is  considered  in  this  section, 

as  a  special  case  of  a  closed  cylinder  (Eq.  (4-55) ).  In  order  to  make  use  of 

.  .  .  .  ...  (22)(24)(33) 

the  extensive  literature  on  diffraction  by  a  circular  cylinder  ,  an 
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elliptic  cylinder  has  been  chosen  whose  principal  axes  are  parallel  and  per¬ 
pendicular  to  the  optic  axis  respectively,  and  related  to  each  other  as  follows: 


b  =  Je 


a  =  1,  e>0 


(4-86) 


where  a  and  b  are  the  half  axes  of  the  ellipse  in  the  x  and  y  directions, 
respectively.  Such  an  elliptic  cylinder  transforms  into  a  circular  cylinder  of 
radius  a  in  the  u-v  space  if  we  let  x=u,  y=Ve  v  (Fig.  20).  While  the  con¬ 
figuration  in  Fig.  20  is  convenient  for  the  relatively  simple  confirmation  of 
the  assumptions  and  results  of  Section  D  leading  to  Eq.  (4-55),  it  lacks  in 
generality  since  the  obstacle  dimensions  depend  on  the  medium  parameter  e. 
This  modd  therefore  does  not  permit  the  independent  investigation  of  the  solu¬ 
tion  as  c  takes  on  complex  values,  as  was  possible  in  previous  sections. 
However,  this  model  is  employed  to  test  in  detail  another  aspect  of  the  pro¬ 
cedure:  the  ability  to  apply  the  scale  transformation  to  the  asymptotic  expres¬ 
sions  in  the  isotropic  region,  and  to  recover  therefrom  the  asymptotic  result 
for  the  anisotropic  case.  The  boundary  value  problem  to  be  solved  is  the  same 
as  in  Eqs.  (4-61)  with  S  given  by  the  equation 


=  1  . 


('  -87) 


Transformed  into  the  u-v  space,  the  boundary  value  problem  to  be  solved  is 


+-^Uk2)G(R,  P.')  =  -6(u-u')  6(v-v')  , 

Sv  / 


ac  .  2  2 

— -  ~0  on  u  +  v  =  a  , 

a  n 

r~z  2 

Radiation  condition  at  R  =  lu  +v  —  “  . 


(4-88a) 


(4-88b) 


(4-88c) 


G  is  defined  by  Eq.  (4-65d),  and  k  =  k(  The  rigorous  solution  to  Eqs.  (4-88) 

is  well  known,  and  various  representations  are  available  which  have  different 
convergence  properties.  The  representation  best  suited  to  our  purpose  is  given 
by<”> 
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Fig.  20  Diffraction  by  an  elliptic  cylinder: 

The  xy  and  uv  spaces  for  observation  point  in  the  shadow 
region  (a,  b)  and  in  the  illuminated  region  (c,  d). 
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G(R,R')  =  jGJp.p'iuJG^tp.cp'juJpdp 


(4-89) 
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(4-89a) 
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(4-89  b) 
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(4-89e) 
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G  (r.c'iu)  =  V  (T2iu)_1  +2TTi)I;  lmn>0  (4-89f) 

v  (j 


hJ|2)(2)  *  n  h(u)(7-) 


(4-89g) 


Eq.  (4-89)  may  alternatively  be  written  in  terms  of  a  series  of  residues  arising 
from  the  poles  of  the  integrand,  which  are  located  at  the  points  »  where 


,(1) 

Hj  (ka)  =  0  .  I 

P 

{ 1 6 )( 2?\ 

These  points  have  been  calculated  by  several  investigators'  ’  and  for 

ka  »  1,  they  are  given  asymptotically  by 


(4-90a) 


U  «kat  •  (ka)  ""  , 

P  P 


(4-90b) 


where  *  are  defined  by  Eq.  (4- 56c),  and  for  the  Neumann  type  boundary 
conditions 
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(4-90c) 


In  view  of  Eq.  (4-89f),  G{p^,^/)  can  be  interpreted  as  the  field  ;.'hich  results 
from  an  infinite  number  of  line  sources  located  in  an  angular  space  of  infinite 
extent  (-®<cp<ao)  and  separated  from  one  another  by  a  distance  of  2TT.  Evi¬ 
dently,  only  one  line  source  lies  in  the  physical  region  0<cp<2n,  and  this 
source  alone  accounts  for  the  geometric  optical  parts  of  the  total  field.  The 
sources  outside  the  physical  region  of  the  angular  space  contribute  only  to  the 
diffraction  effects. 

Eq.  (4-89)  may  now  be  transformed  back  to  the  x-y  space,  (where  it 
is  a  solution  of  Eqs.  (4-6l)  and  (4-87),  and  evaluated  asymptotically.  The 
details  of  this  calculation  are  given  in  appendix  I.  It  should  be  noted  that  the 
solution  as  given  by  Eq.  (4-89)  has  different  convergence  properties,  according 
to  whether  the  observation  point  is  located  in  a  geometrically  illuminated,  or 
in  a  shadow  region.  Thus,  the  asymptotic  evaluation  and  physical  interpretation 
of  the  results  are  different  from  one  another  in  those  regions.  If  the  observa¬ 
tion  point  lies  in  the  shadow  region,  then  the  total  field  is  equal  to  the  diffracted 

(33) 

field.  In  that  case  it  can  be  shown  that  Eq.  (4-89)  becomes 
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(4-91) 


whose  asymptotic  evaluation  for  ka»  1  (after  transforming  back  to  the  x-y 
space)  is  carried  out  in  appendix  I  and  yields: 
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(4-92) 


where  Sj  and  are  the  arc  lengths  and  ^*8*  20a),  respec¬ 

tively,  and  A(q^)  is  the  Airy  function  of  argument  defined  in  Eq.  (4-56d). 
Eq.  (4-92)  should  be  compared  with  Eq.  (4-55). 

To  show  that  Eqs.  (4-92)  and  (4-55)  are  identical,  one  notes  first  from 
Eq.  (4-57a) 
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Thus,  from  Eqs.  (4-90),  (4-92)  and  (4-93) 
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(where  b(-t)  is  the  radius  of  curvature  of  the  cylinder  at  t) ,  and  the  equation 
of  the  cylinder, 

y  =  infeifa^-x^  (4-96) 

it  is  not  difficult  to  verify  via  Eq.  (4-57a)  that 
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(4-97) 
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thereby  confirming  the  structure  of  the  exponent  in  Eq.  (4-55)  with  Eq.  (4-59). 
Moreover,  from  Eqs.  (4-96)  and  (4-57a),  it  is  noted  that 


N  d-t  =  J7  a 


dx 


2  2 
a  -x 


so  that 


(j>N  dl  =  2TTa-/T 


(4-98a) 


(4-98b) 


where  the  (j)  is  taken  around  the  circumference  of  the  cylinder.  Thus,  using 
Eqs.  (4-59),  (4-90b)  and  (4-98), 


i  2 Tin  -  (j)  -^i  kQfN  lv(l) ]  -  a(-t)  j-  dl 


(4-99) 


thereby  making  the  denominators  in  Eqs.  (4-55)  and  (4-92)  equal.  In  a  similar 
manner,  the  amplitude  factor  in  Eq.  (4-92)  may  be  identified  as  the  diffraction 
coefficient  It  has  therefore  been  verified  in  this  special  case  that  the 

asymptotic  evaluation  of  the  scaled  exact  solution  agrees  with  the  result  obtained 
when  the  scale  transformation  is  applied  directly  to  the  asymptotic  formula  in 
the  isotropic  region. 

If  the  observation  point  lies  in  the  illuminated  region,  the  field  consists 
of  a  geometric  optical  part  due  to  the  l  -  0  term  in  Eq.  (4-89)  and  a  diffracted 

part  due  to  the  l  i-  0  terms.  Thus  one  may  write  in  this  region 
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where  ^  stands  for  a  summation  over  all  integers  l  excluding  -6  =  0,  and 

iu|cpV| 

(cp.cp  ;U)  =  - ,2uI-  (4- 100a) 

o 

The  series  part  of  Eq.  (4-100)  is  treated  exactly  as  in  Eq.  (4-91)  and  leads  to 
Eq.  (4-92),  except  that  the  quantities  d^.dj,  d2«  Sj,  S2,  etc.  have  to  be  taken 
now  from  Figs.  (20c)  and  (20d).  The  integral  in  Eq.  (4-100)  furnishes  the 
geometric-optical  field. 

It  has  therefore  been  verified  that  Eq.  (4-92),  which  was  obtained  by 
performing  the  asymptotic  evaluation  in  the  anisotropic  region,  agrees  with  the 
previously  derived  Eq.  (4-55)  wherein  the  asymptotic  result  was  obtained  in  the 
isotropic  (scaled)  system,  with  the  scale  transformation  applied  subsequently. 
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CHAPTER  V 

CONCLUDING  REMARKS 

The  object  of  this  study  has  been  to  gain  an  understanding  of  electro¬ 
magnetic  radiation  and  diffraction  in  anisotropic  dielectric  media.  Emphasis 
has  been  placed  not  only  on  formal  mathematical  solutions  but  on  approximate 
methods  for  the  evaluation  of  explicit  results  as  well  as  on  a  physical  inter¬ 
pretation  of  relevant  phenomena. 

To  achieve  these  aims,  a  special  class  of  problems  has  been  investi¬ 
gated  in  detail:  two-dimensional  problems  of  diffraction  by  variously  shaped 
objects  in  a  uniaxially  anisotropic  dielectric.  By  choosing  the  optic  axis  in 
the  medium  at  right  angles  with  respect  to  the  axis  of  the  cylindrical 
scatterers,  one  finds  that  the  electromagnetic  fields  may  be  determined  from  a 
single  scalar  function,  thereby  reducing  substantially  the  mathematical  com¬ 
plexity  and  facilitating  a  thorough  study  of  the  resulting  formal  solutions.  Two 
parameter  ranges  have  received  special  consideration,  corresponding,  res¬ 
pectively,  to  "small"  and  "large"  obstacle  dimensions.  In  the  first,  the  scat¬ 
tering  properties  are  described  essentially  by  multipoles  of  appropriate 
strength  and  orientation,  while  in  the  second,  the  radiation  characteristics  are 
specified  conveniently  in  ray-optical  terms.  In  each  category,  various  repre¬ 
sentative  problems,  for  which  exact  solutions  may  be  constructed,  have  been 
analyzed  in  detail,  and  the  results  have  been  phrased  in  such  a  manner  as  to 
lend  support  to  the  above-mentioned  physical  mechanisms  which  are  operative 
in  establishing  the  radiation  field. 

A  scaling  method^  may  be  employed  to  relate  solutions  for  the  present 
class  of  scattering  problems  in  the  amsotropic  medium  to  equivalent,  known 
solutions  in  an  isotropic  region  when  the  dielectric  tensor  elements  are  posi¬ 
tive  real.  Therefore,  for  this  parameter  range,  the  construction  of  solutions 
is  elementary,  and  the  remaining  task  has  been  an  "invariant"  phrasing  of  the 
results  so  as  to  include  quantities  which  depend  only  and  explicitly  on  the 
medium  parameters,  and  on  such  structural  features  as  the  radius  of  curvature 
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of  the  obstacle  surface  and  its  orientation  with  respect  to  the  optic  axis. 

Further  and  non-trivial  effort  has  been  required  to  extend  the  range  of  validity 
of  the  results  to  complex  values  of  the  tensor  elements,  thereby  providing 
solutions  for  a  broad  range  of  medium  constants,  including  the  "hyperbolic" 
range  where  one  of  the  tensor  elements  is  negative  real  and  infinite  values  of 
refractive  index  may  occur. 

The  results  for  various  large  obstacle  problems  have  been  compared 
with  predictions  made  by  generalizing  Keller's  geometrical  theory  of  diffraction 
to  anisotropic  media.  Since  that  theory  proceeds  on  a  ray-optical  basis  which 
retains  its  validity  under  quite  general  conditions,  its  confirmation  for  several 
special  structures  admitting  of  an  exact  analysis  provides  a  basis  of  support 
similar  to  that  achieved  for  the  isotropic  case.  It  may  therefore  be  concluded 
that  the  present  study  has  provided  a  general  mechanism  for  the  analysis  of 
two-dimensional  diffraction  problems  in  uniaxially  anisotropic  media  having  the 
afore-mentioned  orientation  of  the  optic  axis,  and  that  the  insight  gained  thereby 
may  serve  as  a  basis  for  the  investigation  of  more  complicated  situations. 
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APPENDIX  A 


Derivation  of  the  Green's  function  representation  in  elliptical  coor¬ 
dinates  (Eg,  (3-11)). 


The  boundary  value  problem  stated  in  Eqs.  (3-8)  is  transformed  via 
Eqs.  (3-9a,  b)  into 
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Radiation  condition  at  2hcosh5  — *  ,  (A-lc) 

where  h  is  defined  by  Eq.  (3- 10a)  and  Fig.  4.  The  homogeneous  equation 
(A- la)  can  be  separated,  yielding 
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,2 
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- ^+(\-2h  cos  2ti)u2  =  0  (A- 2b) 

dq 

with  X  denoting  a  separation  constant.  U^(r|)  has  to  be  periodic  in  T)  with  a 

period  of  2”  .  This  is  one  of  the  properties  of  the  Mathieu  functions  of  the  first 

kind  and  integral  order^,  ce  (n;  h^),  se  (n;h^),  and  me  (r|;h^)  which  is 

mm  m 

related  to  the  former  via  Eqs.  (3-18a,b).  We  note  the  completeness  relation 
of  these  Mathieu  functions, 
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We  may  assume  an  eigenfunction  expansion 
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Substituting  Eq.  (A-4)  in  Eq.  (A- la)  and  using  Eqs.  (A-2b)  and  (A-3)  we  get 
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°m  me  )  m 

m 


(A-5) 


(A-5a) 


The  boundary  conditions  on  g  (?  ,  -3  *)  are 

m 


=  0  at 


^  Jo  ’ 


(A- 5b) 


Radiation  condition  at 


(A-5c) 


Solutions  of  the  homogeneous  equation  (A- 5)  are  the  "radial”  Mathieu  functions* 

M^($;h).  In  particular,  for  our  implied  time  dependence  of  e  *u,fc  ,  the  pair 
m  (1)  {3) 

of  suitable  functions  is  M‘  '(?)  and  M*  '(?),  which  has  the  asymptotic  behavior 

mm 

for  ?-*  , 


M{1)(§;h)-J  (Zhcosh?)  =  J  (fcR)  , 


(A-6a) 


M(3*{§;  h)~H(1){2hcosh?)=  H^*(kR) 
mo  o 


(A- 6b) 


Their  Wronskian  is  given  by 


M^>(S;h>]  =  A  . 


(A-6c) 


If  we  define 


f ,(?)  =  M(3)(?>:h)  , 

1  m  > 


(A- 7a) 


f  (?)  =  m^}{?  ;h) — Mm>(5<;  h)  ’ 
2  m  <  ;h)  m  < 

m  o 


(A-7b) 
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it  is  easy  to  show  that 


W 


2i 

tt 


where  5^,  and  5^  are  defined  by  Eqs.  (3- 11a,  b). 

The  solution  of  Eq.  (A-5)  subject  to  conditions  (A-5b,  c)  is 

(291 

from  and  f  in  the  usual  way  : 


g  (5. 5') 

m 


(A-7c) 

synthesized 

(A-8) 


thus,  from  Eqs.  (A-4),  (A-5a)  and  (A-8),  we  get  Eq.  (3-11). 
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APPENDIX  B 


The  quadrupole  radiation  in  the  field  scattered  by  a  narrow  strip. 
From  Eqs.  (3-19),  { 3 - 20b)  and  (3-20d)  we  see  that 


g^cp.cp')  =  -  sin  2<p  sin  2cp'+ O(h^)  ; 

also,  using  Eqs.  (3-22),  we  have 

TI  .  .  .  -j sin  6  +  -jfcJoi  cos  6  _ _ 

„  _  U  _  u  cos  6+  v  sin  5  _  Ja 1  o  Vc  o  fTT 

scp  =  -B-  = - R -  =  - ^■■MTa-rwrgT -  Vac  = 


(B- 1) 


cosc 


r  N{6o)  N(9) 


e  sin 9  sin0  +  acosS  cos0 
o  o 


N(3o)N{8) 


(B-2) 


_  sin{0  -0)(e  sin3  sin8  +  a  cos  9  cos  8  ) 

sin2qp  =  2  sincp  coscp  =  zjae  - - - - - s - 2-  . 

lN(6)N(3  )]2 

°  (B-3) 

An  expression  similar  to  (B-3),  with  9'  instead  of  S,  is  obtained  for  sin2cp'  . 
Using  Eqs.  (3-17a),  ( B - 1 )  and  (B-3)  we  get 


*2 


“  a  c  ..  . 

sr  (koa| 


^  sin(0o-9)sin(So-6,) 
[N(8)  N(9')]2 


(c  s  in  0s  in  0  +  a  cos  9  cos  9  )  • 
'  o  o 


*  (e  sin  9  ‘  sin  8  +a 
o 


cos 


0' 


cos  c  )  +  0(k  a) 
o  o 


(B-4) 


To  verify  that  Eq.  (B-4)  represents  the  radiation  from  an  electric 
quadrupole  line  source,  and  to  find  the  quadrupole  radiation  intensity  in  the 
scattered  field  (Eqs.  (3-15)),  consider  Eqs.  (3-24a,b)  with  the  source  term 

J  =  £(cos9o  ^+sin9Q  ^:)6(r)  .  (B-5) 

Comparison  of  Eqs.  (B-5)  and  (3-24c)  shows  that  the  source  given  by  Eq.  (B-5) 
corresponds  to  a  linear  quadrupole  as  shown  in  Fig.  B-l. 
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Fig.  B-l  A  linear  quadrupole  as  a  superposition  of  two  dipoles. 


Substituting  Eq.  (B-5)  into  Eqs.  (3-24a,b)  yields  the  following  equation  for  H: 

H  =  -B6(r)  (B~6a) 


where 


3x  By 


B  =  p 


sin  9 


,cos  3ola  •  £^j  ^  cos2so-  c  si"2eol  5^7 


{B-6b} 


If  we  consider  H=BG,  then  the  equation  for  G  reduces  to  Eq.  (3-  8  ),  whose 
solution,  subject  to  a  radiation  condition  at  infinity,  has  already  been  given  by 
G^(£,  o)  in  Eqs.  (3-13).  Thus,  the  solution  of  Eq.  (B-6a)  is 


H(r)  =  |  BH^1  [kQVay2+  C  * 


(B-7a) 


For  k_)rN(S)>>!  one  may  write 
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c  o 


The  derivation  of  Eq.  (B-7b)  from  Eq.  (B-7a)  is  analogous  to  the  derivation  of 

sin(c  -6) 


Eq.  (3-29),  s  _c; 

H(  r)  -  p  a  e  —  — —  (c  sin  3  sin  jq+  a  cos  8  cos  9q)  CkQ  rN(6)]} 


Comparison  of  Eqs.  (B-7b)  and  (B-4)  yields 

1 4 


(B-7b) 


P  =  -A, 


.  „  (k  a)  sin(G  -o'He  sin  5  ,  sin  6  +a  cos  9*  cos  6  )  , 

i£  -°J. - L° - - 2 - —  +0(k  a)6 

Id  i  2  /%  o 


N^G  ) 


J 


(3-8) 


where  A  is  the  intensity  of  the  incident  field  at  the  location  of  the  obstacle. 

Eq.  (B-7b)  shows  that  the  leading  term  of  g2  in  Eq.  (B-4)  is  a  linear  quad- 

rupole  term,  and  Eq.  (B-8)  gives  the  intensity  of  the  quadrupole  induced  in  the 

strip  by  the  incident  field.  It  is  easy  to  see  from  Eq.  (B-7b)  that  the  axes  of 

zero  radiation  in  the  quadrupole  radiation  pattern  are  no  longer  perpendicular 

to  each  other.  While  one  zero  occurs  at  5  =  8  ,  the  other  zero  occurs  at 

tan  9  =  -  —  cot  9  ,  which  reduces  to  5  =  3  --=■  onlv  when  ^  =  1  (isotropic 

C  o  „  Oc  t 

case),  or  when  5  =0.  j  (special  orientation  of  the  strip). 

O  W 
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APPENDIX  C 


Reflection  by  a  perfectly  conducting  plane  in  a  homogeneous  uniaxial 
medium. 

To  derive  Eq.  (4-10)  from  Eq.  (4-1),  consider  Fig.  8.  A  ray  passes 
through  the  point  P(0,  0),  is  reflected  at  the  point  C(x,  y)  and  then  passes 
through  the  point  Q(d,  0).  Then, 

C  Q  C  _  U _ 

J Nd-t  +  jNdf,  =  J l/cos291+e  sinZQl  +  j£Z<\+  e 

P  C  P  1  c  2 


=  Vy2+ex2  +  Vy2+  e(x-d)  =  L(x.  y) 


(C-l) 


The  variables  x,  y  are  subject  to  the  constraint  of  lying  on  the  con¬ 
ducting  plane,  i.  e.  , 

f(x,  y)  =  Bx  +  Ay  -  AB  =  0  (C-2) 

The  equation  which  gives  the  extremum  of  (C-i)  with  the  constrarint  (  C-2)  is 


dL 

y^ 


df 

y^ 


dL 

yy 


df 

*7 


=  A 


•Substituting 


x  .  a 
—  =  tan  h 

y 


x  -  d 


B 


1  ’  y 

one  gets  from  Eq.  (C-3), 

2 


Cx 

_ gjx-.d) _ 1 

J  I7Z  2 
■Vy  +€x 

<-d)2| 

-J 

_.y  _  4 

y 

= 

Vy2*  £x2  - 

/  2  2 
(fy  vtx 

m 

A  ,  a 
•=•  =  tan  y 

B  o 

and  S.  = 
i 

Gj  +  tt 

1  —  V  - 

r  2 


(C-3) 


(tanS.-tane  )  f(e  tan  3  -  l)(tan3. +tan3  )  A  2  tanS  (CtanQ.  tanG  +1)1  = 
i  r  |_  °  i  r  o  l  r  J 


0  , 

( C-4) 
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which,  upon  exclusion  of  the  case  tan  c .  =  tan  9  ,  yields  Eq.  (4-10).  In  this 

derivation  no  restriction  has  been  imposed  on  e  ,  so  that  the  result  is  valid 
also  for  negative  e  .  While  real  rays  propagate  in  this  instance  only  when 


|  tan y 


(C-5) 


Eq.  (4-10)  yields  a  real  6  for  any  3.  . 

r  "  1 

This  is  in  accord  with  an  image  technique'  '*  'by  which  the  reflected 
field  can  be  constructed  as  shown  in  Fig.  C-l.  A  ray  emanating  from  a  source 


S  gives  rise  to  a  reflected  ray  which  appears  to  originate  at  an  image  source 
I.  If  C<  0  ,  incident  rays  can  propagate  only  within  the  wedge  ASB. 


Fig.  C-l  Reflection  by  a  perfectly  conducting  plane. 
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The  image  source  is  then  located  at  the  head  of  the  inverted  wedge  AIB,  where 
the  triangles  ASB  and  AIB  are  congruent  (or  ASBI  is  a  parallelogram,  with 
the  conducting  plane  as  one  of  its  diagonals).  By  using  the  sine  law  in  triangle 
ASB,  one  gets  _ 


SB  _  SA 

sin  (n  -  c  ^"5  )  ~  sin  (8  ^0  j 
co  o  c 


(C-6) 


and  since  AR  is  common  to  triangles  ASR  and  AIR, 


AR 


_  sin(9  -S^) 

^  sin  (tt  -6  +  9  ) 
o  c. 


AS 


sin  (8  +  9  , ) 
c  l 

sin  (9  -  8 , ) 
o  1 


(C-7) 


Upon  combining  Eqs.  (C-7)  and  (C-6),  one  may  derive  an  expression  connecting 
9 1 ,  G^,  0  ,  and  3^  .  If  3^  is  eliminated  by  using  Eq.  (C-5),  one  gets  precisely 
Eq.  (4-10)  relating  3,  to  0  ,  9  and  £.  This  confirms  that  the  reflection 
law  derived  by  Fermat's  principle  agrees  with  the  one  obtained  from  a  solution 
of  the  boundary  value  problem. 

The  field  above  the  conducting  surface  is  given  exactly  by  the  sum  of  the 

contributions  from  the  real  and  image  sources,  regardless  of  whether  8>0 
( t  c 

or  9  1  c  8c  .If  3  j  >  ,  Eq.  (4-10)  still  yields  a  real  8^  .  but  then  neces¬ 

sarily  |  G2I>  |  9cl  ,  because  an  evanescent  incident  ray  giving  rise  to  a  propa¬ 
gating  reflected  ray  would  violate  the  principle  of  conservation  of  energy.  If 
the  incident  ray  coincides  with  SB,  i.e.  then  the  reflected  ray 

coincides  with  IB,  i.e.  0  = -3  .  All  of  this  information  can  be  deduced  either 

r  c 

from  the  analytic  reflection  formula  in  Eq.  (4-10)  or  from  the  geometric  image 
construction.  In  dealing  with  the  conducting  half-plane  problem  (sec.  C)  rather 
than  with  an  infinite  conducting  plane,  the  geometric  optical  direct  and  reflected 
parts  of  the  field  can  be  constructed  in  the  same  way.  If  we  assume  in  Fig.  C-l 
that  the  half-plane^  coincides  with  that  portion  of  the  infinite  plane  which  lies 
to  the  left  of  the  point  R  (i.  c.  ,  the  edge  of  the  half  plane  is  at  R),  then  the 
reflected  ray  RQ  now  constitutes  the  boundary  between  the  region  where  both 
direct  and  reflected  rays  exist  (to  the  ieft  of  RQ)  and  the  region  where  only  the 
direct  rays  exist  (to  the  right  of  RQ).  As  has  been  shown  in  Eqs.  (4-44),  the 
location  of  this  boundary  may  also  be  deduced  from  the  rigorous  solution  of  the 
problem.  The  solution  for  diffraction  of  a  cylindrical  wave  by  a  half-plane  shows 
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that  the  boundary  between  said  regions  has  to  be  real  even  when  the  edge  of  the 
half  plane  is  in  the  shadow  of  the  source  (point  R'  in  Fig.  C-l,  for  example). 
Eq.  (4-10)  confirms  that  this  is  so,  and  that  the  pole  contributions  to  the  fields 
given  by  Eqs.  (4-36)  and  (4-45)  (see  appendix  E)  yield  exactly  the  geometric 
optical  parts  of  the  field. 


6 
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APPENDIX  D 


By  solving  the  boundary  value  problem  of  reflection  of  a  cylindrical 
wave  from  a  perfectly  conducting  plane  (b-“)  in  a  uniaxially  anisotropic 
medium,  Felsen^"^  showed  that  the  reflected  field  is  given  asymptotically  by 


I  H  (P)|  = 

r 


c'4-  d 


(D-l) 


with  0^,  0.,,  d  and  d'  defined  in  Fig.  D-l,  and  Hq  defined  by  Eq.  (4-16). 


y 


Fig.  D-l  Reflection  by  a  perfectly  conducting  plane. 

(The  choice  of  source  Q  and  observation  point  P  in  Fig.  (D-l)  in  a  way  such 
that  -d*  cos  6  ^  =d  cosG^  ,  facilitates  some  of  the  following  calculations  with¬ 
out  loss  of  generality). 

We  want  to  show  that  Eq.  (4-17),  which  has  been  derived  by  employing 
geometric  optical  principles,  reduces  to  Eq.  (D-l),  which  has  been  derived 
from  a  rigorous  solution.  As  in  Appendix  C,  we  write  the  equation  of  the  plane 
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£(x,  y)  =  x-y  tan  6  =0 


(D~2a) 


and  the  optical  path  from  Q  to  P 


L{x,  y)  =  ^(B-y)^+  )*+C  (A-x)^ 


(D-2b) 


It  is  also  shown  there  that  Fermat's  principle  implies 


«,9r6z)  =  0 


Using  the  relations 


S 1  ri  • 


A  - 


1 


yS2+x2 


(D-2c) 


(D-3a) 


cos 


^{B-y)2+  x2 


(D-3b) 


sin 


(A-x) 


*  r 


■y(B-y)^+  ( A-x)^ 


(D-3c) 


cos  v  ,  = 


B-\ 


V^-y)2t  (A-x)2 
one  can  show  from  Eq.  (D-2c)  that 


g(G,,S,)  =  - 3 

6  1  2  cos  o 


N<’V 


NtSj) 


(D-3d) 


Also 


cos  r,1-’  sin  ~  sin  *■  , 
o  2  o  2 


z  sin(r  -$.) 
^  g  _ _ o  1 

^1  N3(9j) 


coT5  cos  9  ,  +  e  sin  9  sin  ST 
o  1  o  1 


=  0 
(D-4) 

(D-5a) 
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e  sin(6  -  9  ) 
_ c  o 

N3(62> 


Thus,  for  the  case  b-“  in  Eq.  (4-17) 

d9. 


dT) 


1 


!  0  =  const 
o 


sin(0o-8  p  N  (02) 

sin(e2-0o)  N^Gj) 


From  Eqs.  (4-14)  one  finds  (for  general  b) 


d9„  d,  sin(0o-02) 


d  S  ^  a  sin(0  ^  -  9  Q) 

(see  Fig.  9).  But  it  is  shown  in  Ref.  (15)  that  for  b-** 

d'N(01)  =  a  N(02)  ; 

thus,  from  Eqs.  (D-5c),  (D-6)  and  (D-7) 


N(62) 


sin(82-90) 


and 


“  sinTB^l 

Imz)d6l  _  /Wep 

/>lm37vzz t=vnto 


(D-5b) 


(D-5c) 


W^2 

for  0  =  constant  or  b-“.  From  Eqs.  (4-12),  (4- 13)  and  Fig.  9  one  can 
o 

directly 


I  H  (P)  |  = 


H 

o 

1 

f  NfS^dS., 

(D-6) 

(D-7) 

(D-8) 

(D-9) 

deduce 

(D- 10) 


which,  upon  using  Eqs.  (D-7)  and  (D-9),  reduces  to  Eq.  (D-l), 
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APPENDIX  E 


Rigorous  solution  of  the  half-plane  diffraction  problem. 

The  method  by  which  1(a)  in  Eq.  (4-35)  is  found,  is  described  n  detail 
in  Ref.  (4)  chapters  I,  II  for  a  half-plane  in  an  isotropic  medium.  Though  the 
following  derivation  is  in  principle  identical  with  the  one  in  the  above  mentioned 
reference,  it  differs  in  enough  details  to  justify  its  presentation  here. 

Consider  Eqs.  (4-22)  for  the  case  of  plane  wave  excitation, 


L(2)cp  =  0 

(E-la) 

=  0  on  Y  =  0;  -»<X<  0 

(E  -  lb) 

with  cp  defined  by  Eq.  (4-25b).  Since  -4>.  in  Eq.  (4-25b)  satisfies 

o 

ii 

__rr 

'NJ 

(E  -2) 

the  equation  for  cpg  become 

=  0 
s 

{ E-3a) 

L“-s-  -L<V,ko(A3VAA)«P[-iko 

b2x]  = 

i  k  €  sin(9  -5') 

= - NTs7) - exp  [-lkoB2xJ  on 

j  Y  =  0 

1  -®<X<  0 

(E-3b) 

Radiation  condition  at  r~® 

(E-3c) 

Edge  condition  at  r- 0 

(E-3d) 

$  (a)  and  i  (a)  defined  by  Eqs.  (4-28)  are  regular  functions  of  a  in  the  upper 

and  lower  halves  of  the  complex  a  plane,  respectively.  (The  analytic  properties 
of  functions  defined  by  integrals  are  discussed  in  detail  in  detail  in  Ref.  (4), 

Chap.  I).  In  subsequent  discussion,  all  functions  with  a  subscript  +  (plus)  or 
-  (minus)  will  be  regular  in  the  upper  or  lower  halves  of  the  complex  a  plane. 
Also,  we  will  denote  b>  .’(X,0  +  ),  l(a,0+)  the  value  of  cp(X,  Y)  or  $(a,  Y)  on 
the  upper  side  of  the  half  plane  (and  its  extension  at  X>0),  and  by  ^(X.O-), 
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$(a,  0-)  the  same  values  on  the  lower  side  of  the  half-plane  (and  its  extension). 


cpg  is  continuous  on  the  surface  Y  =0,  0<X<®;  therefore 

*a*(a.0+)  =  *s+{a,0-)  =  $s+{0)  . 


(E-4a) 


The  tangential  component  of  E_g  is  continuous  on  the  whole  plane  Y=0;  there- 


L{1)*g+(a.o+)  =  L(1)tR+{a,0-)  =  L(1)^s+(0) 


(E-4b) 


L(14  (a.  0+)  =  (a,  0-)  =  lSlh  (0) 

S'*  s  —  s  *• 


(E-4c) 


Using  Eqs.  (E-4)  and  (4-30)  one  finds 


*  (a,0+)  =  i  (a,  0)+?  (a,  04)  =  1(a) 

s  s  •  s  • 


(E-5a) 


*  (a.0->  .(a.  0)4?  (a.0-)  =  -1(a) 

o  S  »  S  ” 


L{1)^s+(0)+ L(1)‘s_(0)  =  ( i  a  A3  -  i  A2q  y )  1(a)  =  i^eA^jA  a%  1(a) 
j^(  1 ) x,  (0)  is  a  known  function*  :  via  Eq.  (E-3b) 


(E-5b) 


(E-5c) 


L(1)^  (0)  =  -jL  j  dx 


r  ■  n  i  n  k  c  sm  8  -9  ) 

f  ,  iax  T  (1).,  .  o  'o' 

dx  e  L'  cp  (x,  0)  =  - == - 

L  s  N(6')^(a-k0B2) 


(E-6) 


Noting  that 


o-aYA2  Vko+  a/N(5oWko'  a/N(9o> 


(E-7) 


We  get  from  Eqs.  (E-5)  and  (E-6) 

L(1)J  (0) 

s+ 


k  €  sin  (9  -9  ') 
o  o 


(ko4a/N(6o)  N(3  ') JBj ko+a/N(eo)(a-koB2) 


_ 


=  TN<%W(9O)US-(0+)45-(0*)]=K'(a)  (E'8) 


*It  is  assumed  that  for  the  class  of  functions  under  consideration  here,  the 

.  t  ( 1)  jfj  i  a  x 
operators  L  and  dx  e  commute. 
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k  e  sin  (8  -8  ') 
o  Q 


k  6  sin  (8  - 6  / ) 
o  o' 


N(6,)j2nJko+a/N(eo)(a-koB2)  N(0  ')V^/k  +  k0B^/N(60,(a'koB2) 


k  €  sin(8  -8 ') 
o  o 


N(0  ')  JST (a-koB2)  [VV'a/N(0j  '  VkD+ 


=  F_(a)+F+(a) 


(E-9) 


K  (a)  and  F  (a)  are  regular  in  the  lower  half,  while  F+(a)  is  regular  in  the 
upper  half,  of  the  complex  a  plane.  Using  Eq.  (E-9)  in  (E-8)  we  get 


^k+a/mel 


+  F,  (a)  =  K  (a)  -  F  (a) 


(E-10) 


The  position  of  the  branch  points  is  shown  in  Fig.  11.  In  Eq.  (E-10)  the  left¬ 
side  is  regular  in  the  upper  half,  and  the  right-hand  side  in  the  lower  half,  of 
the  complex  a  plane,  having  as  a  common  region  of  regularity  the  real  axis. 
Thus  each  side  can  be  equated  to  a  polynomial  P(a),  the  degree  of  which  must 
be  determined  by  the  edge  condition  as  follows: 

We  may  assume  that  near  the  edge  (i.  e.  for  r«  1)  cp^  can  be  repre¬ 
sented  by  a  series  of  the  form  9 

cp  =  P  )  a(S)rn  ;  (E-li) 

s  L>  n 

n=o 


the  total  energy  E  stored  in  the  scattered  electromagnetic  field  in  a  small 
cylindrical  volume  element  surrounding  the  edge  will  be  given  by  an  expression 
of  the  form* 

2"  R  r  (»v  I2  (io  u  '2 


=  J  J 


R 

r 


8=o  r=o  *- 


3^(8) 


\-rt- 


+  a 


12 


tStP=)+ 

f -5T/  +  a22(e)\TTl/  + 


BOX? 


rdrd 
(E- 12) 


In  order  for  E  to  be  finite  for  any  finite  R,  it  is  seen  that  M  in  Eq.  (E-ll) 
must  be  an  arbitrary  positive  quantity.  From  this  result,  the  asymptotic 


*See  footnote  on  p.  41  for  elaboration. 
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behavior  of  $  and  ?  for  i cl }  —  *  may  be  derived,  by  using  the  identity 


J  x^elClx  dx  =  ±r(u  +1)  e  2  ^  *  a41'1 


(E-13) 


It  follows  from  EnS.  (E-8),  (E-9)  and  (E- 10)  that  for  |al-*,  P(a)=0,  and  thus 


K  (a)  -  F  (a)  =  0 


(E-14) 


One  may  find  1(a)  from  Eos.  (E-5)(E-6)  and  (E-14),  the  result  being  Eq.  (4-35). 
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Asymptotic  evaluation  of  the  integrals  in  Eqs.  (4-36)  and  (4-45). 

Equation  (4-36)  is  the  exact  expression  for  the  scattered  part  of  the 
txcic  cue  to  a  piano  wave  wmcn  tails  upon  a  conducting  tiaif-plane  in  a  uniaxialiy 
anisotropic  medium.  It  is  valid  for  all  complex  z  with  lme>0.  If  the 
incident  wave  is  cylindrical  (i.  e.  excited  by  a  line  source  at  finite  distance 
from  the  edge),  the  appropriate  expression  for  the  scattered  part  of  the  field  is 
given  by  Eq.  (4-45).  Both  integrands  have  a  pole,  which  for  certain  ranges  of 
the  observation  point  contributes  a  residue  term  to  the  solution.  The  other 
term  of  the  solution  arises  from  the  integration  over  the  given  path.  For 
h  r  >:>  1.  i.  e.  for  high  frequencies  and/or  observation  points  far  from  the  edge, 
this  term  can  be  calculated  asymptotically  by  means  of  the  saddle  point  method, 
and  given  in  closed  form.  This  conventional  method  fails  and  must  be  modified 
when  the  saddle  point  lies  near  the  pole. 

Via  the  transformation 


a  =  k  N(c  )  sin  ’• 
o  o 


(F-l) 


the  scattered  part  in  Eq.  (<!-  56)  becomes 


r~ “ 

TV? 


ri 

P 


d  l  cos  :  exp 

J  t  i  kQ  r  N(-)  cos  ( x - X' g )J 

-  sin  i 

.  ,  “il 

s,ni  “nF~T 

(F-2) 


where  the  saddle  point 


tan 


is  defined  by 


-V"YX 


1  -  tan 6  tan  - 

o 


(F-  3) 


f7"  Y  (tan  -  -tan  '  ) 

with  the  x  sign  applying  for  Y’O  respectively.  For  <•:  =  l .  Eq.  (F-3)  becomes 


tan  i  .  =  r  tan  .  (F- 3a) 

s  1  '  =  I 

where  (Fig.  10)  is  the  observation  angle  in  the  X-Y  coordinate  system. 

The  pole  of  the  integrand  in  Eq.  (F-2)  is  located  at 


10? 


sin'D 


B2 

NlO 


(F-4a) 


which  expression  can  be  rearranged  by  using  Eq.  (4-29a)  to  read 


tanuu 


1  +  e  tan  6  tan  9  ' 
_ o _ 

P  ^F(tan 8 _-tan S  ') 


(F-4b) 


By  equating  Eqs.(F-4b)  and  (F-3),  one  finds  that  when  a1  =ui  , 

p  s 


tan  9  = 


-(1  -  €  tan  9  )  tan  9  '+  2  tan  9 


( 1  -  €  tan  0  )  +  2€  tan  S  tan  9 
o  o 


.  for  Y>  0 


(F-5a) 


tan 6  =  tan 6  '  for  Y<0 


(F-5b) 


which  limiting  angles  have  already  been  noted  in  Eqs.  (4-44).  The  path  of 
integration  and  the  locations  of  and  Xg  in  the  complex  it)  plane  are  shown 
in  Fig.  F-l(a),  for  €  real  and  positive.  For  other  values  of  t.  the  mapping 
from  the  a-plane  to  the  x-plane  will  be  different.  For  example,  if  e<0  and 


tan  9  > 
o 


VT 


the  picture  in  the  x-plane  will  be  given  by  Fig.  F-l(b).  For  C<0  and 

tan9  < — —  , 

°  Vul 


the  picture  again  looks  different,  but  as  the  associated  calculations  are  analo- 

(15) 

gous  to  the  former,  this  case  will  not  be  discussed  separately. 

It  is  evident  that  if  the  original  path  p  is  deformed  into  the  steepest 

descent  path  SDP  through  the  saddle  point  x  ,  the  pole  at  x  will  be  inter- 

s  p 

cepted  by  the  steepest  descent  path  whenever  ;  this  corresponds  to 

observation  points  lying  in  region  A  of  Fig.  12  if  Y>0,  wherein  a  geometrically 

reflected  part  of  the  field  exists.  When  x  =x  ,  then  the  observation  point  lies 

P  s 

on  the  boundary  between  regions  A  and  B.  The  angle  of  this  boundary  is  given 
by  Eq.  (F-5a)  and  is  recognized  as  the  angle  of  the  limiting  ray  which  is 
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(a) 


(t>) 


Fig.  F-l 

The  original  and  the  steepest  descent  paths  of  integration  for  e  >  0  (a) 

and  e  <  0  (b) 


reflected  geometrically  "at"  the  edge  ■.■{  the  half-plane.  The  pole  contribution 
to  CPg  is  found  by  calculating  the  residue  of  Eq.  (F-2)  at  jy  =  u)  , 


'V  ~  ±  exp 


[±i  k  r  N(")  cos  lx  - 
o  p 


'■S>J 


tor 


Y<  0 


(F-6) 


Upon  substitution  of  the  correct  .L'g  for  Y<0  ,  respectively,  it  can  be 
shown  that  Eq.  (F-6)  gives  exactly  the  geometrically  reflected  plane  wave  in 
region  A  and  that  it  just  cancels  the  incident  wave  in  region  C  cf  Fig.  12.  For 
calculating  the  saddle  point  contribution  to  the  field,  we  use  the  formula 


J 

P 


a  / 

g\ 


'dx 


-2  n 


g  ('« 


f(a,s) 


(F-7a 
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where  in  our  case  (Eq.  F-2) 

:  =  k  rN(9)  (F-7b) 

This  yields  the  expressions  g'  .  en  by  Eqs.  (4-43). 

Eq.  (4-45)  can  be  derived  from  Eq.  (4-36)  as  shown  in  Appendix  G.  The 
pole  contribution  to  Eq.  (4-45)  in  the  region  Y <0  w:ll  be  calculated  as  an 
illustr?'  it.  In  that  region,  Cg  in  Eq.  (4-45)  is 


.  *4  exD  T-i?  X' +  i  a.  (?)  Y  1“  ext>  f-ia  Xx  i  q„(a)  y] 

=  —JL-  dS  - J  I  da  — 1 - ‘± - J 

s  sn2  c  .4  4koW^  -*  ^k0N(60)-a  (a +  5) 


(F-8) 


The  pole  occurs  at  a  -  -§,  and  by  the  residue  method,  the  pole  contribution  to 
the  second  integral  in  Eq.  (F-8)  is  found  to  be 


Upon  noting  that 

tribution  to  -i'  , 
s 

'0 

S  pole 


expfi?X4  i  q ? ( - > ) y! 

Z  T1  i  - L - — - i 

+  5 

q^( -S >  =  -qj(~>)  (see  Eq.  (4-32),  one  obtains  for  the 


-i 

r  exP  j 

d§  -J 

i  ?  ( X-  X  ’ )  -  i  q  x 

<?)(Y-Y'j| 

4-vr 

J  5 

-ao 

Vkow2(8 

i 

®  exp 

|  d? 

fi5(X-X')  -  iqx 

(§)(Y-Y')] 

‘  4" 

j 

-oc 

A2^ql*q2' 

=  H(1)  [k  N(-)  !r-r'jj 

4  o  o  —  — 


(F-9a) 
pole  con- 


(F-9) 


i.  e.  ,  that  part  of  *  which  cancels  the  incident  cylindrical  wave  in  the  geo¬ 
metrical  shadow  sore  (for  verification  of  the  last  step  in  Eq,  (F-9),  see  Eq.  (16) 
in  Ref.  15;  attention  should  be  called  to  the  different  definition  of  ^  anc* 

Aj  A ^  A^  in  that  reference). 

The  asymptotic  evaluation  of  Eq.  (4-43)  by  means  of  the  saddle  point  tech¬ 
nique  is  very  similar  to  the  evaluation  of  Eq.  (4-36),  because  the  equations 


no 


which  define  the  saddle  points  in  the  ro/mplex  a  and  §  planes  are  independent 
of  one  another.  Defining 


0.  =  k  N(9  )  sintu 
o  o 


(F-  10a} 


?  =  k  N(8  )  sin  z 
o  o 


(F-lOb) 


and  substituting  Eqs.  (F-10)  into  Eq.  (F-8), 

j  p  d  z  cos  z  exp  ^i  kQ  r '  N(0)  cos  { z- zs)J 


c c  - 

s 


8r'2JT  J 


p(z) 


<1 1  -  sin  z 


x 


1 


,.  d  v  cos  t  exp  fxi  k  r  N('  )  cos  (u>- uu  ) 
L  o  _ s. 


p(x) 


^1  *  sinu;  (sinu:  t  sin  z) 


(F-  1 1) 


u;  is  defined  by  Eq.  (F-3),  and  z  is  given  analogously  by 
s  s 


a3y'-a?x 

tan  -  -  — — - - — 

•  F  y' 


1  +  f.  tan  9  tan  6  ' 
o 

F<  tan  S  tancQ) 


{F- 12) 


which  is  identical  with  Eq.  (F-4b).  The  simple  saddle  point  integration  breaks 
down  when 

-a*  ~  z  ,  (F- 13) 

s  s 


corresponding  to  observation  points  near  the  shadow  boundary  (see  discussion 
of  Eqs.  (F-3).  (F-4)  and  (F~5)  ).  For  all  other  observation  points,  Eq.  (F-7a) 
can  be  used  separately  for  the  two  integrals  in  Eq.  (F-ll),  yielding  the  expres¬ 
sion  given  in  Eq.  (4-46). 


Ill 


APPENDIX  G 


Derivation  of  Eg.  (4-45)  from  Eg.  (4-36). 

The  synthesis  of  the  line  source  result  (Eq.  (4-45))  from  the  plane  wave 
result  (Eq.  (4-36))  may  be  carried  out  as  follows:  If  the  source  term  in  Eqs. 
(2-2)  is  of  the  form 


M  =  M(5)  =  z  6(Y-Y  ) 


/.  e 


i?  X 


— o 


the  differential  equation  for  cp(?)  becomes 

,2 


A1  “2  +  A2  ^  ZA3  T^Y  +  k'W>  =  ^UW)  =-6(Y-Y')  ^ 


(Aj,  A^,  A^  are  defined  by  Eqs.  (4-24)).  One  assumes  for  <p(?), 

i?  X 


»{X.  Y;  l)  =  h(Y)  e 


which  implies  that 

,2 


d~ 


\d  Y 


A?  dY 


A2/Ti 


(G-l) 


(G-2) 


(G-3) 


(G-4) 


As  v(?)  is  continuous  across  the  plane  Y  =  Y/.  h(Y)  must  aiso  be,  and  from 
Eq.  (G-4)  we  get  by  integration  between  the  limits  Y,-’,  and  Y'+q 


lim 
q  -  0 


dh 
d  Y 


A 


1  2i? 


Y#  +  ■ 


Y' 


-1 


A  i  2 


(G-  o) 


The  solution  to  Eq.  (G-4)  is  constructed  from  the  homogeneous  solutions,  which 

.  -  ,  °*> 
are  given  for  Y  >  Y  by 

io,  A-SHY-Y’)  -iq,  ,U)(Y-Y') 


h(Y)  =  Ce 


LU 


C  e 


5  2.1 


(G-6) 


with  ^(^)  defined  by  Eq.  (4-32),  and  C  being  a  constant.  Imposition  of  con- 
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dition  (G-5)  on  Eq.  (G-6)  leads  to  the  solution  ox  Eq.  (G-2), 

-tq^  (Y-YVi§X 
i  e _ | _ 

2„  A2(q2-1,l) 


<P(5) 


i  e 


iq2,lY'’  -ic^2,lY+UX 


2  V2  n  ^  k  qN(  0O )  TJJkJi(8jn 
It  can  be  shown  that  (see  Eqs.  (4-26) ) 


(G-7) 


exp  j^i?X-iqiYj=  exp  [-ikoB2X- ikQB3Y  ] 


if  we  define 


3 

■3 


-k  B , 
o  2 


(G-8) 


(G-9) 


Then  f»  r  Y'  >  0,  a  plane  wave  falling  onto  the  half  plane  at  Y  =  0  is  given  by 
Eq.  (G-d),  which  is  proportional  to  C0(?)  in  this  region.  The  scattered  field 
cp  (?)  due  to  such  a  plane  -wave  is  given  by  Eq.  (4-36),  which,  u^on  use  of 

sp 

Eqs.  (G-9),  reads 


%(?)  =  * 


VW**  r  X+iqj .  2(a)  y]  [r_  ^ 

^  !  — jrTrasrriia^)  •  lG'10) 


2  n 


J 

-33 


This  solution  remains  valid  for  -®  <?  <®  so  that  the  scattered  field  due  to  the 
source  in  Eq.  (G-l)  is  given  for  Y'  >  0  by 


da  exp^-iaX+iq^  ^(a)  Yj 

) 


iq,(?)Y' 

.  1  rt  UU.  VAl/l  iW  ^  '  ISI,  ' 

*  (?)  =  ±  — is=— . -  ■  f  —  t  ,  ‘»-j 

s  cJkoN(So).?  _i  ^a-koN(0o)  (a +5 


By  applying  the  inverse  Fourier  transform 

1 


(G- 11) 


fix')  =  -~=  j’  d?  e ’ *  *’  X  F (0 ) 


(G-  1 2) 
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to  Eq.  (G-2),  we  get 


(■ 


“T  +  a2^t4  2A3  yj^Y  +  k><x-Y=  X'-Y'> 


which  is  the  desired  form  for  line  source  excitation, 
will  therefore  be  the  inverse  Fourier  transform  of  cp 
is  given  by  Eq.  (4-45). 


=  - 6( X- X ')  6(  Y -  Y  *)  (G- 12) 

The  scattered  part  of  cp 
g(§)  (Eq.  (G-ll)),  which 
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APPENDIX  H 


Asymptotic  evaluation  of  Eq.  (4-80)  in  the  illuminated  region. 


The  Weber  functions  D  j  (t  l  J-2  ik)  and  D  ^  ^  (± r)  ^-2  i  k ) 

'2  +2k  'I+2k 

which  appear  in  the  integral  representation  of  the  field  (Eq.  (4-80)),  satisfy  the 
equations  T  2  1 

-^T  +  k2(C2+>7k2)  D  .  .x(±%<F2ik)  =  0  (H-la) 


‘  2 

~Y  +  k2  (§2+  X/k2)  D  j  u(i§^2lk)  = 

-d  3  J  “  2  +  2k 

'  2  ?  *1 

-2-5-+  k  (r;  ->./k  )  D  ,  .  .  (=  r.  V^iTk )  = 

-dT1  J  "  2  "  2k 


(H-  lb) 


For  high  frequencies  (large  k)  one  may  construct  asymptotic  forms  of  these 
functions  by  the  WKBJ  method  as  follows 


i  .a'5 


*  -1  * 

"  2  2k 


[§2+\/k2] 


exp  i  k  J  -Jx2+  X/k2  dx  (H-2a) 


°  I 

'  2  "  2k 


[-2  i  k)  ~ 


fT - “^174  exP  ik  jVx2-X/k2  dx 

[t2-X/k2J  I  JVk 


(K-2b) 


Di  iX  {'r-PTk)  ~  ; - - ~i74exP  ik  jVx2-x/k2  dx  + 

'2 '2k  -  A / k  j 


- - ^tT7T  exp  -i  k  J  V^-Vk2  dx 

L’*2-  A/kj  -/T/k 


(H-2c) 


2-2 

Assuming  that  §,r  are  real,  Eq.  (H-2a)  is  valid  for  A/k  >  -%  ,  while  Eqs. 

(H-2b,  c)  are  valid  for  > /k2  r  2  .  Outside  these  regions  the  exponentials  become 
real.  For  |A/k'i  »r,2,f  Eqs.  (H-2)  reduce  to  Eq.  (4-75a),  and  the  exponents 
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become  proportional  to  %  and  respectively. 

The  coefficients  A,  B,  C  and  D  in  Eqs.  (H-2)  may  be  evaluated  by  com¬ 
paring  Eqs.  (H-2)  in  the  region  where  p2,  |2  »  |  X/k2{  with  the  asymptotic  forma 
given  by  Eqs.  (4-72). 

It  is  useful  to  note  that 


(H-3a) 


If  we  substitute  the  asymptotic  forms  (H-2a-c)  into  Eq.  (4-80)  we  get 


r  _  if  d* 

i  4k  J  ~  ju,  rr 

--  sln(2k'I> 


-  D  i <?  -/^rrsjD  j  u(-cv^27k) 

1  *'»***  i  —  "t* ’ 


2  2k 


2  2k 


(H-  3b) 


'  x  D  ,  ^  ("  i  k)  D  j  ~ 

~2  ‘  2k  "  2  "  2k 


(H-4b) 
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Eqs.  (H-4)  are  really  valid  only  in  the  region 


W2  r/2 

•mints 


<  <  min 

k 6  l 


{H-5} 


where  min  |a,  b|  stands  for  "the  smaller  of  a  and  b".  Tlius,  if  there  is  a 
real  saddle  point  in  that  interval  one  may  obtain  a  closed  asymptotic  form  for 
Gj  from  Eq.  (H-4a),  which  will  represent  a  propagating  wave.  If  no  saddle 
point  is  found  in  the  interval  (H-5),  the  asymptotic  form  in  Eq.  (H-4a)  becomes 
useless.  Eq.  (4-75a)  indicates  that  in  this  case  the  integral  has  no  real  saddle 
points  at  all,  and  does  not  represent  a  propagating  wave.  Next  it  will  be  shown 
that  the  first  part  of  Eq.  (h-4a)  has  one  real  saddle  point  in  ths  region  (H-5), 
whereas  the  second  part  has  none.  Thus  Gj  is  given  asymptotically  by  the 
first  integral  alone. 

One  may  easily  verify  that 


_d 

dX 


x+yx  ± 

J  2.  X  _  X  , 

c^x  X  t  -J  In — — - 

k  k  /X/ k 


r“r. 


x  +>/x^±- 


r 


=  i  In 


(H-6) 


Thus 


dqj(X) 

sr~ 


=  In 


?w?2+  -4 

k 


5'l  ’■ 


7 


■^7 


(H-7a) 


and  the  saddle  point  condition  for  the  first  part  of  G^  become; 


_s 

Z 


k4 


=  1 


(H-7b) 


This  condition  has  a  simple  geometrical  interpretation  as  shown  in  Fig.  H-l. 
From  Eqs.  (4-66),  it  follows  that 

"  -  J R  +  v  =  VR(  1  +  coscp)  =  ^2R  cos~ 


j  P.  *+  v '  =  Vr'(1*co.p')S  ^R'ccsy- 


(H-8a) 

(H-8b) 
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?W?2+—J  . _  _ _ .  Y-a  .  Y+g 

_ k_  _  Vl  -  costp  4-  vl  -  cosff>  +  2  sin y  sing  _  sin  £  sin  2 

J,  X  +  cos?  +v^l  +  cos<p  -  2siny  sing  cos  4-  cos 

;2._J  22 


=  tan 


and  similarly. 


l;r+Js  +-f 

“  k  Y 

/ - -  T 

'  j.  'Z  S 


(H-10a) 


(H-lOb) 


Vr'  -7 


which  satisfies  the  saddle  point  condition  Eq.(H  -  ?b).  Now  it  is  seen  from  Eq, 

(H-4b)  that  r  • - ^ ~  | - j- -  I - -  1 

ql(Xs}  =  ?yi2+-|+  |§/|^/2+-|  +  T1Jn2--f +HW2— f  I  (H-  11a) 

k  k  *  k  *  k 


- 

*  k  k 


=  Vm  -  cos  cp)>/R(l  -  cosq?)  t  2t  +  VR{1  +  cos<$  V^Rtl+  cos^P)  -2 1  = 

=  R^y/l  -  cos(Y-g)  >/l  -  cos(Y+g)  +  +  cos(v-  a)  ^ 1  +  cos(y+  g)j  =  2R  cosa 


and  similarly 


i'|'Y§*2+_®-  +  ~yr>'2--7  "  ZR'cosg7 


(H-llb) 


(H-llc) 


^us'  "1  P  1  r  “i 

exp  l^jr  q  j  (X  s)J=  expj^ik  (R  cosg  +  R'  cosg')J  =  exp  j^ik  I  _R  -  _R #  M  (H-12) 

which  is  exactly  the  phase  term  at  the  point  P  of  a  wave  emanating  at  Q  (Fig. 
H-l).  The  saddle  point  condition  for  q^X)  in  Eq.  (H-4b)  would  be 


?  4.  Mi  s  \=  1  '  <  I -•  2 ,  ss  •  ,  I  t  2  s 
a  +-y  3  +  — ^  U  1  +J^  + — ^  rt  n  — j 

l'  k  k 


T|  +Jr - y 

¥  k 


(H-13) 


tv|-< 
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Eqg.  (H-8a-d)  show  that  a  condition  like  Eq.  (H-9),  i.  e.  ;  a  X  which  is 
independent  of  <P  and  cp  ,  cannot  be  found  here. 

The  second  part  of  Eq.  (4-80)  yields 


D\  ii<-vC2Tk> 

"  2"  2k 


_L  f  dX  2  2k _ D 

4ki  -44 


D  x  (?v^2ik)  x 


x  D,  iX  {-§  V-2ik)  D  J  a  (rfV-2i  k)  ~ 


*  -L  * 

*2  +  2k 


2  2k 


s{a2bcI  sin(|:|)n[(!2t7)<r2+7,(n''7)!r'a'7i 4 


+  ABD 


where 


1 

4*-  -j—  f(52+4)(?'2t4)(n2-4)(V2-^)Tf  * 

(2k-2lnl  k  k  k  J 


q,M  =  9X(^)  -  2  J  dx  = 
y;./k  1 


XP  -y  q^tM 


(H-  14a) 


J  2  X 

ITT  x  V-Jn°‘7 
=  q,X*)  -  2  nJno- ^  +  2fc2ln  ^ 


(K-14b) 


q x (X)  is  given  by  Eq.  (H-4b). 

It  is  seen  that  the  first  part  of  G2  cancels  (asymptotically)  the  first 
part  of  Gr  It  must  be  noted,  though,  that  Eq.  (H-14a)  is  valid  only  in  the  region 

-min(?2,?'2)  <X  <  r~  (H-15) 

k^ 

This  implies  that  the  geometrical  interpretation  of  Fig.  H-l  and  Eq.  (H-9)  for 
the  first  part  of  G2  is  valid  only  when  the  observation  point  lies  in  the  shadow 
region,  where  t  <  -JZ  .  In  the  illuminated  region  (t  >  no/2),  the  interpretation 
still  holds  for  Gj.  but  the  first  part  of  G2  has  no  saddle  point  when  X/k  >qQ 


and  is  therefore  of  a  smaller  order  of  magnitude  than  . 

The  asymptotic  evaluation  of  the  second  part  of  is  done  as  follows: 
From  Eqs.  (H-14)  and  (H-6)  it  is  seen  that  the  saddle  point  aondition  is 

given  by 


5+J?2+4  rujra  1 


/?,2+ 


s 


X 

s 

-7 
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It  can  be  shown  that  Eq.  (H-16)  is  satisfied  if  X  is  given  by  Eq.  (H-9).  QCP 

s 

is  the  path  of  a  ray  which  emanates  at  Q  and  reaches  P  after  being  reflected 
at  the  surface  according  to  the  laws  of  geometrical  optics.  OS  and  (5B  are 
perpendicular  lines  to  Q'C  and  Cir*  respectively  whose  projections  on  the  v 
axis  are  equal  and  denoted  by  t  in  Eq.  (H-9)  and  Fig.  H-2.  To  show  that  these 
projections  are  vqual,  one  has  to  utilize  the  fact  that  the  reflecting  surface  is  a 
parabolic  cylinder  which  implies  that 

e=I(q>-n)  (H-17a) 

2  o 

and 

a  =  y #;  a#  -  Y  (H-17b) 

o  o 


(0  is  the  angle  between  the  tangent  to  the  parabola  at  C  and  the  v  axis). 
Using  Eqs.  (H-8a-d),  (H-17)  and  the  additional  relations 

cp_ 

=  */R  -  v  =  <JR  (1  -  cosCp  )  =  JZ  R  sin 
o 


=  >/R  -  v  =  */R  (1-  cosCp  )  =  /Z  R  sin  ~ 
o  o  o  o  o  2 

(H-Ha) 

=  ( 

x'+y'  +  tt 

(H-18b) 

=  Y-OL 

(H-18c) 

=  a  +y  +tt  =  a7  +  y'  +** 
o  o 

(H-18d) 

R  sina  sinY 

(H-18e) 

=  < 

R7  sina7  sinY7 

(H-18f) 

R  sina  sinY  -  R  sina7  sinY* 

Loo  o  o 

(H-18g) 

Eqs.  (H-16)  and  (H-9)  can  be  shown  to  be  compatible  by  a  calculation  similar  to 
that  in  Eqs.  (H-10).  A  calculation  similar  to  that  is  in  Eqs.  (H-Tl)  shows  also 

expj^jr  q^(A =  expjik  (R  cos  a  +  R 7cosa'  +  Rq  cosa^-R^  cos  a/)1  =  exp  jik (QC+CP)J 

(H-19) 

which  is  the  phase  of  a  wave  emanating  from  Q  and  reaching  P  by  being 
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reflected  from  the  parabolic  cylinder. 

From  Fig.  H-2  it  can  be  seen  that 

tion  point  P  nears  the  shadow  boundary, 
second  part  of  is  valid  only  in  the  illuminated  region.  Thus  it  is  seen  that 
the  exact  solution,  Eq.  (4-80),  yields  correctly  the  terms  expected  from  geo¬ 
metrical  optics  as  a  first  order  asymptotic  result  at  high  frequencies.  It  would 
be  satisfying  if  we  could  show  that  the  asymptotic  evaluation  of  the  exact  solu¬ 
tion  in  the  illuminated  region  yields  also  the  amplitude  factor  which  has  been 

(24) 

predicted  by  geometrical  optics  (Eq.  (4-17)  .  Jones  claims  that  it  can 
indeed  be  shown  for  the  case  of  a  plane  wave  diffracted  by  a  conducting  parabolic 
cylinder  in  an  isotropic  region,  but  no  detailed  calculation  is  given  in  this 
reference.  At  the  time  of  writing,  this  correspondence  has  not  been  obtained  in 
the  present  study  although  there  is  reason  to  believe  that  the  calculation  can  be 
performed. 

When  e  is  complex  or  negative  real,  the  parameters  §,§,,'n.'n/  and  T] 

become  complex  according  to  Eqs.  (4-81).  If  these  expressions  are  substituted 

into  Eqs.  (H-7b)  and  (H-16),  one  obtains  complex  expressions  of  Xg  as  functions 

of  the  complex  parameter  e.  The  resulting  equations  for  the  loci  of  X  in  the 

s 

complex  X- plane,  as  functions  of  e,  are  extremely  complicated,  and  their 
detailed  study  (for  some  fixed  r  any  £#)  is  probably  done  best  by  computing 
machine- 


2 

t--£- 


or 


X 

s 

7 


t)  when  the  observa- 


The  geometrical  interpretation  of  the 


4 


APPENDIX  I 


Asymptotic  ?■  :;on  of  Equation  (4-89). 

According  to  whether  the  observation  point  lies  in  f;hc  illuminated  or  the 
shadow  region,  there  will  be  two  slightly  different  treatments  '-■?  Eq.  (4-89).  In 
the  shadow  region,  we  follow  Figs.  20(a)  and  20  (b),  assigning  coordinates  to 
the  following  points: 

Q  e  (x'.y'),  Qj  =  (x|,y')«  P  =  (x,y),  P-j  s(x1,y1),  O  =  (0,0)  , 

2  2  2  2  2  2 

(5  3  (u'.V7),  3  (u'.v'),  f  3  (U,  v),  Pj  3  (u^Vj),  "5  3  (0,0)  , 

2  2  2  2  2  2 

The  distance  between  two  points  A  and  B  will  be  denoted  by  (AB);  thus 


p  =  (V5)  = 


=^(u-Ul)2+  (v-v^+u^v2 


?2>  +  (p2! 


f(u-u2)2+(v-v2)2+  tl2+V2 


(1-1) 


p'  =  (53)  =1 


[<5S>)2+  (^(5)' 


/(u#-u  1)“+(v,-v|)+u|  +vj 

^U'U2^+  "  v2^+  U2Z+V22 


d-2) 


a  =  (^Q^  =  (^5Q2)  =  (5^)  =  (C5P2)  = 


_  1  C.  C  1C  C  1  1C  tC  liC  ,1 

-VU1+V1  =VU2+v2  =  VU1  +V1  =  Vu2  +v2 


d-3) 


2n-  lco-cp '  1  -  S  -y  =  2 tt  -|cp-cp'|  -  cos'1  -  -  cos'1 

P  P 

g 

|  cp  -cp '  |  -  3  -  y  =  jcp-cp'l  -  cos'1—  -cos  1  —  =  — 

P  pa 


d-4) 


d-5) 


By  using  the  relations 
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x  =  u  ,  y  = 


F  ' 


one  finds  from  Eqs.  (1-1),  {1-2),  and  (1-3), 


kp/=  2>  +  yi/22+e xi  2 

k  a  --  e  *f  =  =  kjffTT?. 


-  lie,  /£. 

=  k  jyz  +  C  x2 


In  the  shadow  region,  from  Eq.  {4-91) 


r  iU  |c?-(q>  +  2 ”-l)|  _/  7  r  ViU  Icp-tcp'+Zn^)!  _ 

Z  e  P  'U  ^  P 

!.=  -*  \C=  - 1  4=0/ 


d-6) 


=  kQ^('y-y1  2)S  c  (x-x1(  2)*+  y*  2+ e  x^  2  (I-7a) 


(I-7b) 


(I-7c) 


»2nu  -1  (  i |-i  [2~-  |?-T#|]  iu  Icp-tp'h 


•(I-.  P)  e  P"  '  +  e  P  T'j  = 

i2nu  .!  iu  £  -  iu  Lf^l 

=  (l-e  P)  e  PV  2a  A.  P\  2a  /( 


(1-9) 


with  the  last  step  obtained  by  subtracting  Eq.  (1-4)  from  Eq.  (1-5),  and  sub¬ 
stituting  for  | qp-cp y )  .  Sj  and  are  the  arcs  connecting  to  and 
^2  to  in  Fig.  (20b, I  respectively,  and  may  be  expressed  in  terms  of 
Xj  ^  and  X'  2  (see  Eq.  (4-94)). 

The  asymptotic  form  of  the  Hankei  function. 


H0)  (z)  -  /  2  - 

P  irjZ7 

V  p 


Ntt  .  ,2  -iypl 

*  -UP  -  mP  cos  —\  ■ 


(1-10) 


is  valid  for  source  and  observation  points  which  do  not  lie  on  or  very  near  the 
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cylinder  surface.  Eq.  (4-90b)  is  used  for  M  in  Eq.  ( I- 1 0) 

r 


•n  a  "T 

~  —  +  0(k  a)  «  cosy,  as  k  a  »  1 

K  P  P 


(i-ii) 


a  "T 

■r°;  ~  —4  O(ka)  '  ^  cos  8,  as  k  a  »  1 
P 


By  adding  Eqs.  (1-4)  and  (1-5)  we  get 


d-12) 


.  U  .  U  S  +S 

cos  kp-+cos  kf» ~8+Y  =  n-“W1  (I_13) 

and  with  the  relations  derived  above  (Eqs.  (1-7),  (1-9),  (I- 10),  and  (1-13)), 
iumwO)  \°  iHr>-(<P'+2nf,)|  _ 


•  H  1  ^  (kp)  H(1,u  (ko  ')  ^  elMp 


l=-» 


exp  [ik^y-yp^cjx-x/'t  ik^y'-y jpp  ^-] 
ko[{y'yl)Z+ €{x_xl)2]^  [(y'-yl)2+  c(x'-xj)?']^  [l-e  r] 


+  g<x-x2)  +ikp^y/~y/2)  ~  g(x' -x^)  +Mp  — 

ko[(y'y2,2+  C(x-x/]4  ty/-y2)2f  4  [1_e  *1 


(1-14) 
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The  coefficient  which  appears  in  Eq.  (4-91) 


T~ 


Hj(2)(ka) 


is  identical  with  the  one  appearing  in  the  corresponding  isotropic  problem 
(diffraction  by  a  circular  cylinder  in  an  isotropic  medium),  and  has  been  evalu¬ 
ated  asymptotically  for  large  ka  by  various  authors^ 

For  the  illuminated  region,  we  follow  Figs.  (20c)  and  (20d)  for  various 
definitions  and  note  that  and  are  now  the  arcs  connecting  ^  to 
and  (X,  to  respectively  (Fig.  (20d)).  They  are  interpretable  as  paths  of 

creeping  rays  which  shed  from  the  surface  toward  the  observation  point  P 
located  in  the  illuminated  region.  Thus  the  arc  S 2  is  now  the  long  arc  con¬ 
necting  the  points  Q 2  and  P2  in  Fig.  (20d).  From  the  same  figure  one  notes 
that  Eqs.  (i-4)  and  (1-5)  have  to  be  replaced  by: 


2  TT  -  |  Cp  -Cp  '  | 

-  3  -  Y 

=  2n  -  |  cp  -cp '  | 

- 1  a 
-  cos  — 
P 

-  cos-  1  2  = 
p 

a 

(1-15) 

2n  +  lo-qj7 

|-3-v 

=  2  n  +  \v-y'\ 

-1  a 
-  cos  — 
P 

- 1  a 

-cos  — T  - 
p 

fz 

a 

d-16) 

Proceeding  as  above,  one  may  then  derive  a  result  like  Eq.  (1-14),  except  that 
the  -6=0  term  is  not  included  in  the  summation.  Thus  the  expression  given  by 
Eq.  (4-92)  (or  equivalently  Eq.  4-55)  is  valid  both  in  the  shadow  and  the  illumi¬ 
nated  regions,  but  in  the  illuminated  region  one  has  as  well  the  integral  in  Eq. 
(4-100). 
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